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This paper applies the Boltzmann method of gaseous kinetics to the prob- 
lem of charged particles moving through a gas under the influence of a static, 
uniform electric field. The particle density is assumed to be vanishing low, 
and the ion-atom collisions are assumed elastic, but the field is taken to be 
strong; that is the energy which it imparts to the charges is not assumed 
negligible in comparison to thermal energy. In Part I, the formal framework 
of such a theory is built up; the motion in the field is describable by the drift 
velocity concept, and the smoothing out of density variations as an aniso- 
tropic diffusion process. In Part II, the "high field" case is treated in detail; 
this is the case, for which thermal motion of the gas molecules is negligible; 
the equation is solved completely for the case that the mean free time between 
collisions may be treated as independent of speed; complete solutions are 
also presented for extreme mass ratios of the ions and the molecules; special 
attention is given to the case of equal masses, which has to be handled by 
numerical methods. In Part III, information about the "intermediate 
field" case is collected; with the help of a convolution theorem the case of 
constant mean free time is solved; beyond this, only the case of small ion 
mass (electrons) is available. In Part IV, the diffusion process, whose 
existence was proved in Part I, is pushed through to numerical residts. 
Part V discusses the scope of the results achieved and demonstrates the possi- 
bility of extending them semiquantitatively beyond their original range. 

Part I — General Theory of Strong Field Motion 

IA. QUALITATIVE DISCUSSION 

It is well known that if we consider a mixture of gases under no ex- 
ternal forces the steady velocity distribution which establishes itself 
in the mixture does not depend on the interactions between the gas 
molecules; we have always a Maxwellian distribution for each species 
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with a temperature common to all. This result arises from statistical 
mechanics; the derivation of it is simple and requires few assumptions, 
yet it enjoys a wide degree of generality. As soon, however, as a non- 
equilibrium feature is imposed upon the system this simplicity vanishes, 
and the subject acquires ramifications. Results must now be derived by 
kinetic theory. The amount of labor required increases, while, at the 
same time, the result achieved becomes less general. 

A mixture of charged particles (ions or electrons ; in the following often 
simply referred to as ions) and gas molecules can in principle never be 
in equilibrium since the presence of the former in itself represents an 
instability. However, one might expect, that equilibrium exists in a 
restricted sense, for instance, as regards motion. Even this is rarely the 
case under actual conditions of observation. The non-equilibrium fea- 
tures of greatest importance for analyzing ion motion are a constant 
force (electric field) acting upon one species but not the other (mobility 
theory), and a concentration gradient for one particular species (diffusion 
theory). It is the purpose of this paper to apply kinetic theory to these 
problems, and to compute with its help the most important properties 
which such a gas of charged particles possesses. The work will be dis- 
tinguished from similar ones in that the electric field will not be sup- 
posed weak; velocity distributions which have no resemblance to the 
Maxwellian distribution will thus make their appearance. Furthermore, 
the mass of the charged particles will not be assumed small, which 
means the possibility of getting results for gaseous ions as well as elec- 
trons. Magnetic fields, plasma and A.C. phenomena will, however, be 
excluded. The quantities of interest under those conditions are the drift 
velocity of the ions, their energy, energy partition and diffusion con- 
stants. These quantities will be calculated by assuming plausible mechan- 
ical models. The work just outlined has been published in part in 
abbreviated form in the Physical Review; 1 the exposition to follow will, 
however, proceed independently from these articles. 

Much of the work which concerns itself with transport processes in 
gases makes use of perturbation theory. This method permits us to 
predict the behavior of a gaseous assembly under an electric field or a 
concentration gradient in the limit when the field or the gradient are 
vanishingly small. The result of so perturbing a Maxwellian distribution 
can be expressed through certain constants, such as the mobility or the 
diffusion coefficient, which involve the Maxwellian distribution and 
the internal interactions, but not the perturbation itself. 



1 Wannier, G. H., Phys. Rev., 83, p. 281, 1951 and Phys. Rev., 87, p. 795, 1952. 
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The limits of such a procedure can easily be estimated. In the case 
of an electric field, perturbation techniques apply if the kinetic energy 
acquired by the ion from the field is small compared to thermal energy. 
This means at least that the energy acquired in one mean free path be 
small, i.e., 

eE\ « kT 

where e is the electronic charge, E the electric field, k Boltzmann's con- 
stant, T the absolute temperature, and X the mean free path. Actually 
the situation is not even that favorable. If the mass of the ions and the 
molecules is very different, the energy transferred upon collision is 
small, and hence the ions possess the ability to store the acquired energy 
through many collisions; for this reason, the inequality reads more 
properly 



M . rn\ 
m\ + M) 



. eEX « kT, 
A 

where m is the mass of the ions and^M the mass of the gas molecules. 
After some substitutions this estimate becomes 

K + ™)eE« P ,, (1) 

m M/ 

where p is the true gas pressure and a the collision cross-section. Taking 
as an example an ion travelling in the parent gas we find 

or in commonly employed units 

F 

- « 2 volt/cm (mm Hg) . 

V 
It is clear that this limit is often surpassed in experimental situations. 

The cases in which the limit (1) is applicable are of no further interest 
here because they are well covered in the literature. 2 A field will be called 
"low" when it satisfies the criterion (1) and "high" when the inequality 
is reversed. It is important to notice that a fixed field at a fixed gas 
density may shift from "low" to "high" through a drop in temperature. 

All calculations to follow will contain the assumption of "low ion con- 
centration" which is often made in studies of this sort. It means that 

2 See for instance: A. M. Tyndall, The Mobility of Positive Ions in Gases, 
Cambridge University Press, 1938, Chapter IV. 
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all effects which ions exert upon each other are neglected. The equation 
for the distribution function of ionic velocities is then linear instead of 
quadratic. It is clear that this simplification presents great advantages 
from the point of view of calculation. 

In deriving a criterion for the validity of this assumption we must 
distinguish two types of effects of the ions upon each other. The first 
is the space charge effect. In this effect the ions at large distances make 
the major contribution. Its magnitude depends on apparatus dimensions. 
The criterion for no space charge distortion of the field E is 

4neX 

where n is the number density x of the ions and X a suitable length chosen 
from apparatus dimensions. Inequality (2) is quite stringent because 
it predicts field distortions at values of n of the order of 10 cm" . This 
is the value at which it will become impossible, or at least difficult, to 
make significant experimental measurements. But from the point of view 
of theory this criterion is not relevant. Space charge does not change 
the character of the velocity distribution of the ions because the type 
of ion-ion interaction producing the space charge field is long range and 
creates only a smooth modification of the electric field which we may 
presume to have been included in the original field. What we are con- 
cerned with here are ion-ion interactions which have a random character 
and thus are apt to upset a velocity distribution derived from the "low 
concentration" theory. From this point of view neighboring ions are 
most effective because their relative location fluctuates rapidly, and 
hence, the Coulomb force between them will induce mutual scattering. 
The magnitude of this force is of the order e 2 n 2/3 where n is the number 
density of the ions. It is known from theory 3 that the effect of a Coulomb 
force is preferably not represented by discrete "collisons" but by a 
continuous bending of the entire path. Thus we come to the conclusion 
that random ion-ion forces have no effect if the force given above cannot 
produce a significant deflection in one mean free path. This means 

e 2 n 2/3 X « mean ion energy (3) 

According to whether we are in the high or low field region we get differ- 
ent criteria from this. At low field the thermal energy predominates 
and we get 

cV " « pa (3a) 

3 Mott and Massey, The Theory of Atomic Collision, Oxford Press 1933, 
Chapter III. 
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At high field the "field" energy predominates and we get 

eV" « eE (M. + » ) (3b) 

A rough evaluation of inequality (3a) for one mm Hg pressure gives 

n <K 10 10 particles/cm 3 

This corresponds to a current of about 10 15 particles/cm 2 sec or 
200 (uamps/cm 2 . At lower pressure the criterion becomes more stringent. 
Equation (3b) gives similar results. 

It is appropriate to survey at this point the past theoretical work 
treating the "low concentration" theory of ionic motion for arbitrary 
fields. A rather complete body of work exists for electrons where the 
following three assumptions seem appropriate: (a) that the mass of an 
"ion" is very small compared to the mass of a molecule, (b) that the 
total kinetic energy is conserved in each encounter, and (c) that the 
angular distribution is isotropic in the center of mass system. 

These three assumptions lead to a distribution law given by Chapman 
and Cowling. 4 The law has considerable flexibility because it permits 
the substitution of an arbitrary relationship connecting mean free path 
and speed of encounter. In addition it contains no assumption as to 
whether we have low or high field. A more specialized and explicit dis- 
tribution law is obtained if we assume in addition: (d) that the colli- 
sion cross-section is independent of the speed of encounter (hard sphere 
approximation) ; and (e) that we deal with the high field case only. The 
special law resulting in this case is the distribution law of Druyvesteyn. 

If an improvement over the Chapman-Cowling distribution for elec- 
trons is desired account should be taken of inelastic collisions, that is 
assumption (b) should be discarded. Work in that direction has been 
carried out by Smit, Allen 5 and others. 

The assumption to be discarded first in theory of ionic motion is, of 
course, assumption (a). In order to understand what this implies we 
must understand what advantages assumption (a) has in a calculation. 
In the limit when the ionic mass is very small the encounters with gas 

* Chapman-Cowling, The Mathematical Theory of Non-uniform Gases, Cam- 
bridge University Press 1939, Sections 18.7-18.74. Other references are found 
there. 

« Smit, J. A., Physica, 3, p. 543, 1937 and H. W. Allen, Phys. Rev., 50, p. 707, 

1937. 
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molecules become such that momentum is lost quickly, but energy is 
accumulated in the form of random motion. As a result of this we end up 
with a distribution function which is very nearly spherically symmetrical 
in velocity space. Such a situation permits obvious procedures through 
which the entire calculation is simplified. These procedures will not longer 
be available when assumption (a) is dropped. 

Knowledge concerning the structure of the velocity distribution func- 
tion for gaseous ions is practically nonexistent at this time. Hershey, 
who deals with the motion of ions in the high field case, simply substi- 
tutes for it a Maxwell ian distribution with an unknown offset of the 
origin and unknown temperature parameter, 6 shown in Fig. 1(a). He 
then computes these two parameters by applying the laws of conserva- 
tion of momentum and kinetic energy. It is to be expected that this 
procedure should give reasonable values for the mobility and the mean 
energy of the ion; indeed, if we consider the polarization force only, we 
get exactly the right values; the reason for this is that one may evaluate 
velocity averages for inverse fifth power forces ignoring the distribution 
function and that he did this in effect for the drift velocity and the 





c x c x 

(a) (b) 

Fig. 1 — Simplified pictures for the high field velocity distribution of gaseous 
ions, (a) Hershey's assumption, (b) Modification with correct second moments. 

6 Hershey, A. V., Phys. Rev., 66, p. 916, 1939. 

7 This will be shown in Section IIB. 
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total energy. In order to test whether an offset Maxwellian distribution 
is a satisfactory approximation we have to go one step further and 
examine the partition of the energy among the three degrees of freedom. 
There we find Hershey's distribution in error, for he assumes equiparti- 
tion for the random motion, while, in reality, the random energy parallel 
to the field is much higher than at right angles, 7 giving the distribution a 
decided "ridge" structure. This discrepancy could be taken into account 
by the use of an elliptically distorted Maxwellian distribution, shown in 
Fig. 1(b), and this may prove to be convenient in some applications. 

For a detailed knowledge of the distribution function it is necessary 
to specify the interaction between an ion and a molecule. This interaction 
can be, broadly speaking, summarized under three headings: (a) the 
polarization force, (b) the short distance repulsion, and (c) symmetry 
effects. The polarization force arises because an ion, when passing close 
to a molecule, induces on it a dipole moment; this moment is then at- 
tracted by the charge of the ion. The attractive force F resulting from 
this is 

F = 2 -^ (4) 

where P is the polarizability of a gas molecule and e the charge of the 
ion. The force varies inversely as the fifth power of the distance p; for 
such a force the cross section a varies inversely as the speed of encounter 
7. Whenever the cross section shows this type of variation it is advantage- 
ous to define a mean free time t rather than a mean free path X. The 
formula is 

T = J- (5) 

Nay 

There is a standard difficulty which arises when one tries to make use 
of a formula of the type (5). For most force laws, a total cross-section a 
cannot be defined; a differential cross section per unit solid angle always 
exists, but it becomes infinite in the forward direction because of small 
deflections suffered by particles passing by each other at a large distance. 
Thus equation (5) is, strictly speaking, meaningless. This is actually 
never a difficulty in the computation of a physical quantity. However, 
equation (5) is convenient for order-of-magnitude thinking and the 
question arises how it can be reasonably interpreted. The general method 
of salvaging (5) — excluding a small forward cone from consideration — 
is of little value for this purpose. An analysis of the inverse fourth power 
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attractive potential shows a better way out. The potential gives rise to 
two kinds of orbits; orbits of large angular momentum which look some- 
what like hyperbolas, shown in Fig. 2(a), and orbits of small angular 
momentum for which the particles are "sucked" toward each other in a 
spiralling movement until a repulsive force reverses the trend, as shown in 
Fig. 2(b). A calculation of Hass6 9 shows that the latter type of motion 
is much more efficient in scattering than the former and one gets there- 
for a picture which is semiquantitatively correct if one substitutes into 
(5) the cross-section for spiralling collisions and assumes isotropic scat- 
tering. 10 This cross section equals 



= 2 *1/I+I 



Vp 



(6) 



A numerical estimate of the cross section (6) automatically leads one 
to compare it with the short distance repulsion familiar from the kinetic 
theory of gases. The two are of the same order, but for the usual gaseous 
speeds (which enter into (6) through 7) and small molecules the cross 





ib; 



Fig. 2 — Sample orbits (schematic) showing the motion of a particle in the 
polarization force field, (a) Hyperbolic orbit (large angular momentum), (b) 
Spiralling orbit (small angular momentum). 

8 There are quantum mechanical analogues to these classical ideas; they should 
lead to practically identical answers unless the angular momentum quantum 
number is small. 

9 Hasse, H. R., Phil. Mag., 1, p. 139, 1926. 

10 This will be discussed more fully in Section IIIB. 
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section (6) is bigger. This situation is accentuated in an actual scattering 
calculation which shows an attractive force to be generally more efficient 
than a repulsive force of equal range. 

A detailed numerical discussion of these questions is found in Massey 
and Mohr 11 for the case of He + ions moving through He gas. Their 
interest is in the low field mobility. They show that for this problem the 
repulsive force makes so little difference that it could be neglected en- 
tirely without much affecting the results. It does finally come out that 
the polarization force gives a mobility which is too big by a factor of two. 
But the additional scattering is due to an effect which we listed above 
under (c) : namely a resonance attraction between the He atom and the 
He + ion for which the cross section is abnormally large. It should be 
possible to eliminate this effect by increasing the cross section (6) until 
it masks even this special effect. Lowering the field is not sufficient to 
achieve this because of the temperature motion; it would be necessary 
in addition to reduce the absolute temperature by a sizeable factor and 
so to decrease the value of 7 in (6). Thus we are led to the prediction 
that if the temperature of He is reduced the mobility of He ions in He 
should gradually rise from its "anomalous" value of 12 cm 2 /volt sec to 
the "normal" value of 22 cm 2 /volt sec, which one gets by taking account 
of polarization forces only. 

IB. GLOSSARY 

The complicated appearance of equations in gaseous kinetics suggests 
special care in the use of symbols and a convenient arrangement for the 
reader to find their meaning. It is hoped that the glossary to follow will 
accomplish this purpose. It explains all symbols except those used at 
one location only. 

Generally, Latin capital letters will refer to the gas molecules and 
Latin lower case letters to the ions; Greek letters will have no special 
relationship; exceptions will be made for generally recognized symbols. 
Thus we define 
E, E = electric field. 

x,y = cartesian coordinates at right angles to the field direction. 
z = cartesian coordinate along the field direction, 
r = position vector with components x, y, 2. 
t = time, 
m = ionic mass, 
e = ionic charge. 

11 Massey, H. S. W., C. B. 0. Mohr, Proc. Roy. Soc, 144A, p. 554, 1931. 
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a, a = — = ionic acceleration. 
m 

b = impact parameter. 

bum = limiting value of the impact parameter separating hyperbolic 

and spiralling orbits (equation (125)). 
c, c, c', c a , Cf , C{ = various ionic velocities or components. 
e x , e y , e z = energies of ionic motion (or "high field" parts thereof) 

along x, y, z. 
e* = random part of the above energy, 
j t = total particle current density of the ions (may be a function of 

r and t). 
j = partial current density induced by the concentration gradient; see 

equation (22). 
k = Boltzmann's constant (only when followed directly by T). 
k, k = relative concentration gradient of the ions (a different use is 

made of k in Section HE). 
n = number density of the ions (may be a function of r and t) . 
P> Q> r > s = undetermined constants; used three times independently 

(equations (82), (90) and (160)). 
P ( °\ V = various approximations to these numbers. 
u, u', v = ionic velocities. 

w = ionic velocity rendered dimensionless (see eq. (75) or (85)). 
£f = the inner integral in the double integral eq. (69). 
C, C = molecular velocities. 
2D = ionic diffusion tensor. 
D\\ , D± = components of above tensor parallel and perpendicular to 

the field. 

M = molecular mass. 

iV = number density of the molecules. 

P = molecular polarizability. 

T = gas temperature. 

U, U' = molecular velocities. 

X, Y = left and right hand sides of equation (Ilia). 

1 
/3 = rr^; = temperature parameter (A different use of /3 is made 

in Section IIIB where it is the relative impact parameter b/bu m ). 
Y> y'j It v' = relative velocities of ion and molecule. 
£, rj, f = cartesian coordinates oriented on c. 
p = distance between ion and molecule, 
o- = collision cross section of ions and molecules (may be a function of 7). 
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X = — - = mean free path of ion between collisions with molecules 
No 

(may be a function of 7). 

= mean free time for the ion between collisions with mole- 



~" Nay 

cules. 
t, = same parameter for "spiralling" collisions. 

a = + 1. It is assumed constant in Section ID. 

dlny 

x, Xc Xu = angle of scattering of ion and molecule in the center of 
mass system. 

k = angle of scattering of the ion by a molecule in the laboratory system. 

e = scattering azimuth of ion and molecule in the center of mass system. 

to = scattering azimuth in the laboratory system (azimuth of the initial 
ion velocity about the final ion velocity). 

#, #' = angle between velocity vector and field direction. 

^, ,p, 6,<f>,8 = other angles (these angles are defined on spherical tri- 
angles which are exhibited in Figs. 8 and 15. 

d( c> r> t) = density function of ions in phase space. 

m(c) = ( — ) exp (—f3mc) = Maxwellian velocity distribution function 
for ionic mass. 

^(C) _ {£=. J exp (-/3MC 2 ) = Maxwellian velocity distribution func- 
tion for molecular mass. 

h(c) = "high field" distribution function of the ions for the case that 
the spatial distribution is uniform (the exact meaning of this 
term is to be explained in the text). 

f( c ) = true velocity distribution of the ions for the case that the spatial 
distribution is uniform. 

g(c) = correction to /(c) or h{c) for the case of a constant relative con- 
centration gradient k. 

5(c) = vectorial 5-f unction in velocity space. 

Ei(x) = I — d£ (suppression of two minus signs). 

Jx ? 

I Q (x) = modified Bessel function of order 0. 

K (x), Ki(x) = Modified Hankel functions of order 0, i. (Alteration of 

2 
Macdonald function by a factor — ). 

P,(x) = Legendre Polynomials. 
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h„(c), g y (c) = expansion coefficients which result when h(c), g(c) are 
expanded in Legendre Polynomials about the field direction. 

LAx) = A set of functions of the scattering angle defined in (48). 

( ) = the quantity in pointed brackets is to be averaged. 

(s, v) = abbreviation for {w'P p (cos t?)); the average is taken over /i(w). 

{s, v) = A normalized correction to (s, v) contributed by #(w); see 
equation (155). 
A special convention will be adopted to distinguish velocities before 

and after a collision: 

c', C = velocities before the collision. 

c, C = velocities after the collision. 

When used in this fashion the twelve components of the four vectors 
above satisfy the four identities: 

mc' + MC = wc + MC (7) 

mc' 2 -f- MC' 2 = mc 2 + MC 2 (8) 

The same convention is to apply to other vector quadruples, such as 

u, U, u', U' 

For the velocities in the center of mass system we use 

y' = c' — C = relative velocity before the collision. 

Y = c — C = relative velocity after the collision. 

In consequence of (7) and (8) the 7's obey the relation 

y 1 - t 2 0) 

The multiple integrations occurring in the theory are of the following 
two types. Either they are over the three components of a velocity in 
a Cartesian velocity space; we shall denote such integrations by 
dc, dn, dU', etc. Or they are proper "collision" integrations which classi- 
cally have the form 

7 b db dt 

where b is an impact parameter and e an azimuth. In most cases these 
integrals depend on extraneous factors for their convergence but this 
fact is usually disregarded for convenience; we shall follow this habit by 



• 
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writing the above differential in the form 

— 7o-(7)Il(x) sin xdxde = — 7cr(7)ll(x) dQ y 
4ir 47r 

Here dQ, is meant to represent an integration over a solid angle and the 
subscript 7, that it is over the solid angle swept out by the vector y. 
The notation makes use of the fact that the choice of the polar axis is 
arbitrary in such an integration. The function II(x) is the probability 
of scattering which equals unity for isotropic scattering. In cases where 
small angle scattering is infinitely probable the above expression becomes 
meaningless, strictly speaking, n(x) being a S-i unction at x = and a 
being infinite. However if a quantity such as 1 - cos x is multiplied 
in, which removes the 5-function then the integration gives a finite num- 
ber which may be denoted by (a-(l — cos x))- 

IC. FORMAL SURVEY OF THE THEORY 

Under the assumptions stated in Part IA we may describe the motion 
of ions in a gas by their density in phase space. The change in time of 
this function is described by a Boltzmann equation 12 which, in our 
notation, reads 

dd(c, r, t) dd(c, r, t) c dd(c, r, t) 

dt * ac ' L dr (1Q) 

= ^- ff {M(C')d(c, r, t) - M(C)d(c, r, t)}ya(y)U(x) dOy dC 

The equation is linear in the unknown function d(c, r, t) ; this is due to 
neglect of ion-ion collisions, as stated earlier. The negative term on the 
right hand side actually reduces to a known function of c multiplying 
d(c, r, t). The positive term is a genuine integral term; it has been shown 
by Pidduck 13 that the number of integrations in it can be brought down 
from five to three; this reduction will not be made use of in the following. 
If there were no terms on the left hand side of equation (10) then the 
solution of it would have the equilibrium form 

die, r, t) = nm(c) (11) 

where n is a constant. This result is a direct consequence of equation 
(8) which makes the curly bracket in (10) vanish identically when 
Maxwellian functions are inserted. 

11 See Reference 4. 

13 Pidduck, F. B., Proc. Lond! Math. Soc, 16, p. 89, 1915. 
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The function m(c) is not the solution of our problem because of the 
presence of the second and third term on the left which arise from an 
electric field and a density variation respectively. These disturbances 
will be assumed of different relative importance. The density variation 
will be assumed sufficiently small so that the third term can be treated 
by perturbation theory; the field term, on the other hand will be taken 
so large that the equilibrium distribution (11) no longer represents a 
first approximation to the solution. In consequence, the equation is 
solved in two stages. In the first, only the second term on the left is 
retained, and the resultant equation is treated rigorously; in the second, 
the full equation (10) is used, but the new terms are taken as 
perturbations. 

The first stage describes those properties of the ion gas which it po- 
sesses when assumed of uniform density. Since the field is also assumed 
uniform and not changing in time, the dependence on r and t drops out. 
We may then write 

d(c, r, t) = nf(c) (12) 

where n is a constant and /(c) is a velocity distribution function. The 
equation for / reads 

& 'fc = LII < il/(C ' ) - / * (c ' ) " ^(C)/(c)}7<K7)n(x) da, dC (13) 
with the side condition 

/(c)dc = 1 (14) 



/ 



As a result of solving (13) we shall obtain the distribution function 
/(c) as a function of the electric field contained in a. This distribution 
differs essentially from the Maxwellian one in that it is not symmetric 
about the origin. The vectorial mean of the velocity is therefore not zero 

<c> = //(c)c dc * (15) 

This is the drift velocity of the ion in the field which is reached as a 
compromise between the acceleration a and the frictional losses caused 
by the ion-atom collisions. From the structure of equation (13) there is 
one general prediction that can be made concerning this velocity, namely 
that it depends on the gas density and the field only through a/N; 
this is the well known E/p of the experimental analysis. This type of 
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dependence does not only hold for (c), but for all averages derivable from 
/(c), notably the mean energy. 

A more important formal prediction can be made about the second 
stage of the contemplated calculation. For it will be shown now that the 
diffusion concept is still applicable in the presence of a strong electric 
field. It is true, that if we have a variable density in space the primary 
motion observed is not diffusive but a displacement of the entire density 
pattern with the drift velocity (c). However, once this dominant com- 
ponent is subtracted out, then a supplementary current proportional to 
the density gradient is identified. The constant of proportionality is 
anisotropic, that is, we have a diffusion tensor rather than a diffusion 
coefficient. The tensor is axially symmetric about the field direction, 
yielding a longitudinal and a transverse diffusion coefficient. 

To demonstrate these features it is convenient to assume a special 
type of variation of ion density in space. As we shall see the velocity 
distribution is primarily sensitive to the relative density gradient k; 
we shall therefore assume it to be a constant. In other words we set 

n(r,t) = no exp [k-(r - (c)t)] (16) 

The relation can of course not hold everywhere since n increases beyond 
all bounds in one direction, but we must remember here that we are 
doing perturbation theory, that is k is assumed small. The inconsistencies 
in the assumption (16) can then be pushed as far away as we please. 
Furthermore there is no inconsistency at all in the half space where n 
decreases. It is to be observed that according to the assumption (16) 
the spatial distribution is moving unchanged through space with the 
drift velocity (c). This seems to contradict the program of finding the 
effect of diffusion upon n(r, t). However, we follow in this simply conven- 
tional steady state computational methods in which a gradient is as- 
sumed maintained from an infinitely strong source; the modification 
appears then as a change in the velocity distribution function, which, in 
turn, yields a steady diffusion current. We set therefore 

d(c, r, t) = n(r, t)\f(c) + 0(c)] (17) 

where g(c) is a correction to the solution of (13) which arises from the 
assumption (16). It follows from the definition of n(r, t) and (14) that 

/ g(c)dc = (18) 

The consistency of the assumptions (16) and (17) with equation 
(10) becomes evident when they are substituted into this equation. We 



(19) 
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find, after simplification with (13) 

a "^ + £ // W^ ~ W)0(c')}7*(Y)n(x) d$V dC 

= -k-(c-(c)){/(c)+g(c)} 

This is an equation in velocity space only, r and / having disappeared 
completely; this justifies the assumptions. In solving the equation we 
observe that our interest is only in diffusion, that is, the current resulting 
from a concentration gradient when treated in first order perturbation. 
In this case both k and g(c) are to be treated as small and their product 
in (19) is to be neglected. The equation then becomes 



^*T + S // {M(c)g(c) ~ WMOb^Mx) da,, dc 

= -k.(c-<c»/(c) 



(20) 



The homogeneous prototype of this inhomogeneous equation is (13) ; an 
arbitrary amount of /(c) could thus be added to a particular solution 
of (20) were it not for the orthogonality condition (18) which makes the 
solution definite. 

The existence of the diffusion phenomenon follows easily from equation 
(20). The total current j< is given by 

j,(r, t) = J d(c, r, t)cdc (21) 

Upon substitution of (17) into this expression two terms result 

j«(r, = n(r, 0(c) + j(r, I) (22) 

with 

j(r, = n(r, t) f g(c)c dc (23) 

The first term in (22) is seen from (15) to just equal the product of the 
density and the drift velocity; this is the expected drift current. The 
new current j(r, t) induced by the density gradient is thus given by (23). 
From (20) it follows that ^(c) is a linear function of the three com- 
ponents k x ,k u , k t with coefficients which do not depend on the density 
or its gradient, but only on the unperturbed velocity distribution /(c) ; 
furthermore, the first two of these coefficients are equal. Hence, from 
(23) j comes out as a linear function of the three quantities n(r, t) ■ k x , 
n(r, t)-k v , n(x, t)-k t ; these are the components of the density gradient 
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as is evident from (16); in addition the multipliers of the first two com- 
ponents are equal. We may write therefore 

fc4--C0f£ < 24 > 

where (2>) is a tensor which is axially symmetric about the field direction; 
its two components which we shall call the longitudinal diffusion coeffi- 
cient Dn and the transverse coefficient Dx are computed entirely from 
the unperturbed velocity distribution /(c). This makes D\\ and D± inde- 
pendent of the density or its gradient; it is to be noted, however, that 
they do depend on the electric field as a parameter because this quantity 
enters several times in the course of the computation. 

ID. DIMENSIONAL ANALYSIS 

Dimensional analysis is a convenient tool in a qualitative discussion 
of (13) and (20). In order to get results the situation has to be 
schematized somewhat, but not so much as to impair its usefulness. 
In the first place it is convenient to keep in mind the two limiting cases 
of high and low field, as discussed in the introduction. In addition some 
assumption must be made about a(y) and H(x) occurring under the in- 
tegral sign. The most convenient way to dispose of II (x) is to take it as 
independent of y. This happens to be true for the two models treated 
in detail later, the polarization force model, and the hard sphere model. 
Actually II (x) can be taken as approximately independent of 7 in a wider 
sense. The forces which produce scattering are either repulsive or short 
range attractive, that is, long range attractice forces are absent. As long 
as this is the case the scattering is roughly isotropic and hence can 
change but little with 7. 14 

A more drastic assumption is needed to dispose of a(y). We must 
assume 

a( y ).y a = r (25) 

where a and r are taken to be constants. This assumption contains two 
important special cases in it. They arise respectively by taking a = 
and a = 1. The case a = is the case of a constant mean free path as 
exemplified by the hard sphere model. The case a = 1 is the case of 
constant mean free time; it is applicable to the polarization force as dis- 

14 This statement is checked in detail in Section IIIB for the polarization 
force. This is the attractive force with the longest range which can arise in this 
field. 
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cussed in Section IA. When (25) is inserted into (13) it is seen that 
a, N and r enter only in the combination a/NT. The quantity 

a \_L 

— ]2-« 

NT/ 

has the dimension of a velocity. A second such quantity is 

which arises from the Maxwellian functions under the integral sign. 15 In 
the high field case, this quantity does not enter, that is, the velocity dis- 
tribution functions for the molecules could be replaced by 5-functions 
at the origin. Hence the first combination controls all velocity averages. 
For the mean drift velocity, we can thus write 

(c.) = const '{JftF* (26a) 

This formula gives the variation of the drift velocity with the electric 
field. It is worth while writing the result out explicitly for the two special 
cases discussed above. The first is the case of constant mean free path, 
a = 0, for which 

<c.) = const- a 1 ' V 2 (26b) 

This is a drift velocity varying as the square root of the field or a mo- 
bility varying inversely as the square root of the field. The second case 
is the one of constant mean free time a = 1, for which 

(c t ) = const -ar (26c) 

This means a drift velocity proportional to the field or a constant mo- 
bility. 

In the low field case we cannot disregard one of the two velocity 
parameters constructed above; but now equation (13) is to be solved by 
perturbation theory only, it then yields a drift velocity varying with the 
first power of a/NT. Dimensional analysis then yields the dependence of 
the mobility on the temperature. We find 

(c.) = const-^ (g)T (27a) 



16 Dimensional analysis is incapable of distinguishing between m and M; this 
means that we cannot master dependence on mass by the method of this section; 
all our "pure numbers" are actually unknown functions of m/M. 
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with the special cases 

<C) = const -aX (^) '" (27b) 

for constant mean free path and 

<c,) = const -ar (27c) 

for constant mean free time. Comparison of (26c) and (27c) might lead 
one to surmise that we have here twice the same formula. This is indeed 
the case, as will be shown in Section IIIA. 

Proceeding now to the diffusion problem, we observe from (20), (23) 
and (24) that we must add the quantity 

nk 
NT 

to the previous list of parameters when computing the diffusion current. 
However, the current is always linear in this quantity, which means 
that the diffusion coefficients contain the factor 

J_ 
NT 

and beyond this factor depend on the same variables as previously. This 
gives in the high field case 

D - const.i^p (28a) 

with the special cases 

D = const- a 1/2 X 3/2 (28b) 

and 

D = const -aV 3 (28c) 

In the low field case, the diffusion process becomes independent of the 
field and we get 

D = const ~ (IQH* (29a) 



st 'i\rr \m) 

with the special formulas 

D = const-X^Y 2 (29b) 
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and 



U = COnSt 'T -rr-p 

M 



(29c) 



The information in the formulas (29) is now new, but dependent on (27) 
through a universal relation first discovered by Nernst and derived 
independently for gases by J. J. Thomson; it is widely known as the 
Einstein relation. It states that 



D m d(fi ' )kT 

da 



m 



(30) 



Equation (30) contains of course more than is obtainable from (27) 
and (29), since it relates one undetermined constant to another in a 
known way. 

The dimensional methods of this section are convenient for a rough 
classification of experimental material. Figs. 3 to 7 show the drift velocities 
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Fig. 3 — Drift velocity in an electric field of He + ions in helium gas. Compari- 
son of observed results with an "asymptotic" straight line of slope }$. 
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NUMBER DENSITY 

Fig. 4 — Drift velocity in an electric field of Ne+ ions in neon gas. Comparison 
of observed results with an "asymptotic" straight line of slope \i 

in the parent gas, observed for He+, Ne+, A+ Kr+ and Xe+. The plot 
is a log-log plot of these quantities against a/N, a variety of fields 
having been used to determine each point. The data are taken from 
measurements of J. A. Hornbeck 18, 17 and R. N. Varney. 18 These data 
verify in the first place that the drift velocity depends on a and N only 
in the combination a/N. Beyond this we see that the curves consist of 
two straight line portions: in the lower field portion (c.) is proportional 
to a/N, in the higher is proportional to y/a/N. We recognize in this 
latter region the high field dependence predicted in equation (26b). 
We learn from this that the collision cross section between noble gas 
atoms and their ions is approximately constant in the experimentally 
significant velocity range. To determine these collision cross sections 
the computation of only a single number, namely the one entering into 

» Hornbeck, J. A. and G. H. Wannier., Phys. Rev., 82, p. 458, 1951. 
» Hornbeck, J. A., Phys. Rev., 84, p. 615, 1951. 
18 Varney, R. N., Phys. Rev., 88, p. 362, 1952. 
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Fig. 5 — Drift velocity in an electric field of A + ions in argon gas. Comparison 
of observed results with an "asymptotic" straight line of slope \&- 
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Fig. 6 — Drift velocity in an electric field of Kr + ions in krypton gas. Com- 
parison of observed results with an "asymptotic" straight line of slope >^. 
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(26b) is required. The linear range of this plot is not as informative as 
the high field one. The slope unity is common to all formulas (27), and 
the temperature dependence of the mobility is needed to give the correct 
interpretation with the methods developed here. There is a certain 
likelihood that the parameter a of equation (25) drifts from to 1 as 
the speed of the ions is reduced; this was pointed out for the special case 
of He + in He in section IA. A qualitatively similar situation appears to 
prevail for the other noble gases. 

Part II — The Motion of Uniform Ion Streams in the High 

Field Case 



HA. FORMULATIONS OF THE BOLTZMANN EQUATION 

The dimensional analysis of the last section shows that there is an 
intrinsic simplicity to the high field case which is comparable to the low 
field case, while the intermediate case is more difficult. With one excep- 
tion, however, theoretical analysis has occupied itself with the low 
field case only. We shall try to remedy this in the following. To begin 
with, a tractable but accurate formulation of the problem has to be 
found. Such a formulation cannot treat the field term of equation (13) as 
a perturbation term, but must try instead to make use of the basically 
simple features of the problem, notably those exhibited by the dimen- 
sional analysis of Section ID. 

The equation governing the high field properties of the ions is ob- 
tained simply by substituting 5-functions for the Maxwellian velocity 
distributions in equation (13). This gives 

+ -L A(c) =~ ff «(CW) -^r n <*) d ®y dc ( 3l ) 

OCz t{C) 47T J J T\C) 

A reduction of the number of integrations from five to two must be 
possible in the integral term of (31), owing to the presence of the 5-f unc- 
tion. To achieve this we must transform the variables of integration so 
as to make three of the differentials equal to dC. We do this in the 
following way. First observe that 

Y = c - C 

and that c is a constant vector. Hence we may replace dC by dy. The 
five-fold integration reads then 

dQ y > dC = y dy dti y dQ y > (32) 

that is, it goes over the magnitude y which the two vectors have in 



dh(c) 
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common, and their orientations, for which they are independent. It is 
known that in integrating over the two angles defining an orientation 
the polar axis may be chosen freely. We shall, in the following, adopt c as 
our polar axis with ^, k being pole distances of y and y' to c and <p, w the 
corresponding azimuths. In Fig. 8 these angles are exhibited on the 
unit sphere. All vectors are assumed to be plotted from the center of the 
sphere, and show up through their piercing points. The angles between 
the vectors then show up as sides and the azimuths as angles. The ex- 
pression (32) becomes then 

y 2 dy sin \f/ d\f/ d<p sin k die du 

The main transformation consists now in introducing the three com- 
ponents of C in the place k, $ and tp. The transformation formulas 
follow from the vector identity 



C ' = c — r — T' (33) 



and read in full 



C < " "ITS Sin * C0S V ~whn Sin ' ° 0S " 

, My ■ . • my ■ 

C' = — sin \p sin <p — 57-r — sln K sin w 

' M + m Y M + m 




x 

Fig. 8 — Definition of the angles employed in the formulations of the Boltz- 
mann equation for the high field case. 
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n' My , my 

M + m M + m 

From these equations the value of the Jacobian comes out to be 



mM 2 y 



sin if/ I cos ^ sin k — sin \p cos k cos (ip — oi) } 



(M + m) 3 

We needs its value only at the position C$ = C v = C[ = 0. If we take 
the above equations for C[ , C' n , C' t and multiply them respectively by 
cos k cos to, cos k sin o>, — sin k, add and set C = we get the identity 

My 



M + m 



{ cos >// sin k — sin <p cos k cos (<p — a>) } = c sin k 



The curly bracket is exactly the one occurring in the Jacobian which 
therefore reduces to 

J , ' . % - = — — c c sin ^ sin k 



and hence 



7 2 d 7 dfi r da y , = (M + m ) ^c' dc <h> 
Mmc 



Substituting finally this expression into (32) and (31) we get the Boltz- 
mann equation in the form 

dh(c) . 1 , , v 
a — ^ + -^ h(c) = 

dc z t{C) 

M +m (34) 

= (M + m) 2 f\ M ~ m \ ' 1 



r\M—m\ 1 . r" 

/ -in* n(x) rfc' / h(c') do, 

Jc t(c) Jo 



k-wMmc Jc t(c') 

The equation is in need of additional elucidation as regards the exact 
meaning of c' as a vector and as regards the auxiliary variable x- As to 
the first point we may describe the integration as occurring over a sur- 
face in velocity space. This surface is obtained from the relation 

C = c' = Y ' (35) 

which substituted into (33) becomes 

(M + m)c - my' = My (36) 
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Squaring this and using (9) we get 

(M - m)c'- + 2mc'c - (M + m)c 2 = (37) 

This is the equation of a sphere in velocity space which passes through 
the point c' = c. For all other points c' is bigger than c (collision with 
a stationary object always brings energy loss). The center of the sphere 
lies on the line joining c to the origin; it lies on the side of the origin 
from c when m < M , at infinity (making the sphere a plane) when 
m = M and away from the origin when m \> M. We make use of (37) 
to express the polar angle k of c' with respect to c (which does not occur 
as an integration variable in (34)) in terms of c' '. We get 



cos K = 



(M + m)c 2 - (M - nQc* 
2mcc' 



(38) 



The angle of scattering in the center of mass system also results from 
squaring of (36) if the term my' is first taken to the right. We find 

cos x , W+jnf I _ M> + m> m 

C0SX 2Mm c' 2 2Mm 

There is a more useful form of equation (34) which results if x is taken as 
one of the integration variables rather than c'. Substitution is made 
from the equation (39) above; it yields 

. m + 1 HO = i /" sin x dx ^ (fj f He') A. (40) 
dc t t(c) 4tt Jo t{c ) \c / Jo 

The magnitude of c' and its polar angle with respect to c are now auxiliary 

parameters; the first is obtained from (39) 

M + m 



' C y/W + m 2 + 2Mm cos x 
and the second from (38) and (41) 



(41) 



m + M cos x f*n\ 

C0S K = y/M* + m 2 + 2Mm cos x 
As previously, the azimuth w of c' about c is an independent variable. 

The simplifications of the equation (31) exhibited in (34) and (40) 
still leave a double integral in the fundamental equation. The integration 
over du will now be eliminated by decomposition of /i(c) in spherical 
harmonics about the field direction. There is no loss of generality in 
this step. 

00 

h(c) = X) h,(c)P v (cos #) (43) 
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We have now to consider simultaneously the three vectors c, c' and a 
as well as the angles between them. These angles are denned in Fig 8. 
We study equation (34) or (40) term by term in order to see what be- 
comes of it upon substitution of (43). Starting with h(c') under the inte- 
gral sign we get from Fig. 8 and the addition theorem for spherical har- 
monics 

00 

h(c') = £ K(c')[P, (cos d)P, (cos k) 

+ 2 i: ^T^j P: (cos *)** (cos k) cos mo)] 
„-i [V + n) ! 



For this expression, the integration over <o is elementary and gives 

»2t oo 

/ h(c) da = 2w X) K(c')P, (cos &)P, (cos k) (44) 

Jo »=o 

Further, we get for the derivative in (34) or (40) 

4~ (E h> (fi)P, (cos *)) 

m £ dhAci i_ {(v + 1} (cog &) vp ^ ■ 

,-o «c zv + 1 



+ - KKfi) ^ t ! } {^-i (cos 0) - P^i (cos *)} 

Through the equations (43), (44) and (45), all terms in equation (34) 
or (40) are developed in spherical harmonics with respect to the angle i? 
between c and the field direction. We can therefore annul separately 
the coefficient of each Legendre polynomial in cos #. This gives the follow- 
ing set of equations 

(M + m£ ti£k'hJ*) Pt (cos ,) u(x) # _ hM 
2Mmc Jc t(c) t(c) 

va (dhy-xifi) v — 1 



2v - 1 V dc 



A_i(c) (46) 



+ -27TT V~&- + -7" K+l{c) ) 
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or 

i r htf) /VY Pr (cos k) n(x) sin x dx 

2 Jo t(c') \c / 

_ /Oc) = »o f ^v-i(c) _ Llli k^U)) (47) 

t(c) 2r - 1 \ efc c J 

( y +l)g fjUg) + v+^ fc (c) \ 
+ 2v + 3 \ dc T c ' +U 7 

where 

v = 0, 1, 2, 3 • • • . 

The auxiliary parameters entering are given by (38) and (39) for equa- 
tion (46), and (41) and (42) for equation (47). 

The equations (46) or (47) obtained by Legendre decomposition still 
are, in general, mixed integral-differential equations in one independent 
variable. Further simplification is possible only in special cases some 
of which will be discussed later. An even more simple and tractable form 
of the Boltzmann equation can be achieved in general, however, if one 
gives up the idea of determining the velocity distribution function and 
concentrates instead on its moments. In other words, the Boltzmann 
equation can be looked upon as a system of relations between velocity 
averages, and as such it becomes a linear algebraic system. 

To carry out this reduction we multiply equation (47) by c and 
integrate from to ». The second term on the left is then a simple 
velocity average. The same is true on the right hand side if two integra- 
tions by part are permissible and leave no integrated out part, s ^ — 1 is 
probably adequate for this. The integral over the integral term at first 
looks as follows 

- l f c' +2 dc [' ^ (-Y P, (cos k) n(x) sin x dx 

2 Jo Jo t(c') \c ) 

In this expression we pass from c to c' as the independent variable. 
From (41) we see that 

dc _ dc 
7" V 

Hence the expression becomes 

r c u + 2 wA dc r i r (c\ Pf (cos k) n(x) sin x dx 

Jo t{c ) 2 Jo \c / 
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From (41) and (42) it is seen that this is actually the product of two 
independent integrals if the angular distribution II (x) is independent 
of the velocity of encounter c'. The first integral is then identical with 
the one arising from the second term in (47), and the second is a collision 
integral having no connection with the velocity distribution. Even if 
this is not the case, the second integral is still a dynamic average which 
can be evaluated as a function of c' previously to any knowledge of 
h{c'). We express this by introducing the abbreviation 

I.., = feY P, (cos k) 
Using (41) and (42) we see that /,,„ is the following function of x 
r f v /V-M 2 + m 2 + 2Mm cos xY 

1 8, AX/ — I TJ 1 J • 

\ M + m J 



■■(- 



i# (48a) 

m + M cos x 



V-M 2 4- m 2 + 2Mm cos x 
which, for the particular case of equal masses, takes the simple form 

Z..,(x) - cos' \ x P. (cos |x) (48b) 

With this definition the integrated equation (47) reads 

f / l^jm \ ,. (c)c - dc m ±+ ! + 1) r ^ l(c)c - & 

Jo \ (it{c) I 2v — 1 Jo 



+ 



(y + l)(j - v) 



[ h >+1 (c)c' +1 dc 

Jo 



2» + 3 
or in terms of averages 

(2v + 1) ( 1 ~ 7 /-' (x) c'P, (cos *)) 

= »0> + 8 + lXtT'P^i (cos *)) (49) 

+ ( y + l)(s - w)(<T*Pwn (cos t?)> 

I believe that equation (49) contains all possible derivable relations 
between averages as special cases. Some of the most notable ones are 
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listed below 

s = 1, v = 1 

/l - cos x A Af + m (t . n s 

( — — c cos # ) = — rj — v.o0a; 

\ ar(c) I M 

s = 2, v = 

Z 1 ~ cos * c 2 \ = (M + m)2 (c cos 0) (50b) 

\ ar(c) I Mm 

s = 2, v = 2 

/3Msin 2 x+4m(l-cosx) ^ (cos ^\ = 

X ar(c) ' (50c) 

= 4(MW (ccost?) 

While the averages entering into (49) are not always the desired 
ones, it remains true nevertheless that all solution methods evolved in 
the following use this equation system as a starting point rather than 
other forms of the Boltzmann equation. 

HB. THE MEAN FREE TIME MODEL AT HIGH FIELD 

If the angular distribution in the center of mass system is independent 
of speed and the collision cross section varies inversely as the speed 
then the developments of the previous section permit actually a solution 
of the Boltzmann equation. It is a solution in the sense that all signifi- 
cant velocity averages can be obtained directly without the knowledge 
of the velocity distribution function. 

Before developing these facts from the equations of the last section, I 
should point out that the derivation to follow is in a sense artificial. It 
has been shown already by Maxwell 19 for related problems that if the 
mean free time between collisions is assumed constant specially simple 
techniques may be employed to get constants of experimental im- 
portance. These techniques can be employed here; they consist essentially 
in multiplying (13) by a suitable multiplier, followed by integration 
over c. However, if we were to follow this procedure we would have to 
duplicate for a special model in an unsystematic way the work done 
systematically for all laws of interactions in the preceding section. A 
further advantage of using systematic procedure is that we can see at a 

19 Maxwell, J. C, Collected Papers, Vol. II, p. 40. 
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glance what averages can or cannot be obtained, and what the relation- 
ship is between the high field and the general averages. 

For this reason we limit ourselves at present to the high field averages 
obtainable from (49) . For the special case under discussion this equation 
system takes the form 

(2, + 1) ( 1 ~ J ' M ) (c'P, (cos *)) 
\ at I 



= v{v + s + lXc'-'P^i (cos *)) 

+ (v + l)(s - pXc^P+h (cos *)) 



(51) 



that is we have a system of linear relations connecting the averages 
(c'P, (cos #)). The connection between these averages is made apparent 
in Fig. 9. Each average (c"P v (cos &)) is marked in this figure as a dot in 
an s-v-planc if s is integer. The equations (51) connecting these averages 
are shown as lines with different equations leading to the same dot shown 
in different outline. These equations generally have the shape of a V\ 




Fig. 9 — Interconnection established by the Boltzmann equation among the 
averages (c'P, (cos #)); case of constant mean free time. 
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there are two notable exceptions to this rule, however, which make the 
recurrence method possible, the equations v = have no left leg and the 
equations s = v have no right leg. Starting out with the average 
s = 0, v = 0, which equals unity by definition one can thus proceed 
systematically as shown in Fig. 10, to get other averages. The averages 
reached are the ones for which s and v are non-negative integers of equal 
parity with the restriction s ^ v. One verifies easily that this set is equiva- 
lent to the set of all products of integer powers of the velocity compo- 
nents. 

The first three relations one uses in the path outlined in Fig. 10, are 
the simplified forms of the three equations (50). We find 

2 4(M + m) 3 

(C P 2 (COS &)) = / m giri 2 x + 4m(1 _ CQg x) v i _ CQS g v 

\ ar A O-r I 

or, more conveniently with the help of (53) 

,._ . x , IM sin 2 x + 4m(l - cos x)\ 

(M -f- m) ( / 

(ft = \ 21 I— (54) 

Kt/ .„ /3M sin 2 x + 4m(l - cos x ) \/ l - cos x\ 2 

M 2 m ( A 7. / 

\ ar /\ ar / 

The three equations (52), (53) and (54) give the drift velocity, the 
total energy, and the energy partition of the travelling ion. Equation 
(52) gives a constant mobility and can actually be derived from a low 
field theory. Formula (52) thus states that for problems involving a 
constant mean free time the high field and low field mobilities are numeri- 
cally identical. One would suspect that the intermediate field value would 
have to fall in line too. This is indeed the case as will be shown in Section 
IIIA. 

A convenient interpretation of (53) may be had by combining (52) 
and (53) in the following way 

(mc 2 ) = m{c,f + M{cf (55) 

The left side is essentially the total energy of the ion, the first term on 
the right is the energy visible in the drift motion; it follows therefore 
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that the second term is the "invisible" or random part of the mean 
energy. Formula (55) thus states that 



random energy _ molecular mass 
visible energy ion mass 



(56) 



that is, it exhibits in a quantitative way the capacity of storing energy 
in the form of random motion which light ions travelling in a heavy gas 
possess; for ions travelling in the parent gas the ordered and the random 
part of the energy are just equal; for heavy ions in a light gas the dis- 
ordered fraction becomes negligible. 

There are various ways of understanding the implications of equation 
(54) . One way is to derive the mean energy in a direction at right angles 
to the field by the use of (53). We find 



(d) = 

Mm 



\ ar I 



/3M sin 2 x + 4m(l - cos x)\/ 1 - cos x \ 2 



(IT 



(57) 



\ ar A 

Now from (54) and (57) the partition of the energy e may be obtained. 




Fig. 10 — Order to be followed in computing by recursion the averages 
{c'P, (cos 0)); case of constant mean free time. 
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It comes out to be 

_ / M sin 2 x\ /Mam 2 x\./ Af sin 2 X + 4m(l - cos x) \ / g gN 
e.:e,:*. - ^ - /-\ T /"\ r / 

This shows up immediately the equipartition property for small m/M 
and the overwhelming preponderance of the motion in the field direction 
for large m/M. As an analysis of the random motion, however, equation 
(58) is deficient because the three directions become only comparable 
after the square of the drift velocity (52) is subtracted out of the z-com- 
ponent. We find 

,_. . v 2 /M sin 2 x + 2m(l - cos x) 5 
(Af + m) 



( r 2 \ -(c) 2 - N 91 - (59) 

" »# /3M" sin X + 4m(l - cos x) \/ l - cos x \ 2 

Mm ( X / 

\ ar /\ ar I 

and from this a more refined partition formula which only counts random 
motion 

* /sin 2 x\ /sin 2 x\./M sin 2 x + 2m(l - co s x)\ / fi0 % 

*:*:* = \-7-/ : \-^/'\ JM + «)r / (W) 

For small m/Jlf this result does not differ essentially from (58) but if 
m/M is large the z component of the random energy does not grow in- 
definitely the way the total energy does. Instead it stops at a value which 
is about four times one of the other two values. 

A discussion of these expressions for special models will be delayed 
until the equations are extended to intermediate field conditions. This 
will be done in Part III. 

IIC. THE CASE OF LARGE MASS RATIOS*. ELECTRONS OR HEAVY IONS 

The distribution of velocities for a small value of m/M is treated 
in the literature because it applies to electrons. 4,6 However, for the sake 
of completeness the derivation will be carried out here for the high 
field case. In this derivation all features of the law of scattering are left 
open, except that conservation of the kinetic energy is assumed. 

The development in spherical harmonics carried out in Section IIA 
is the suitable starting point for small m/M, because, in this case, the 
distribution is almost spherically symmetrical and the expansion in 
spherical harmonics is also an expansion in powers of m/M . If we keep 
only ho(c) and hi(c) in the system (47) and treat m/M as small we get 
two equations, one for v = and the one f or v = 1. We may then use 
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(41) and (42) in the simplified form 

<T ~ X + M (1 " cos x) (61) 

COS K ~ cos x + Yf sm x (®^) 

and develop the equations in powers of m/M. Starting out with the 
simpler equation v = 1 we find 

1 - cos x\ , M dfo(c) , , 
5-v — )/ti(c) = -— r - (63) 

This same procedure is not adequate for the equation v = because 
the two left hand terms in (47) cancel in zero approximation. We must 
therefore develop the integral up to linear terms in m/M; this is a per- 
fectly straightforward, though somewhat cumbersome, step. It leads to 
the following equation 



m 
It 






The equation may be integrated by multiplication with c ; this yields 

Wc) = 3^(I^2L*W) (64) 

Elimination of hi(c) from (63) and (64) gives a differential equation for 
Ih(c) which is easily solved by quadrature; the result is 

Except for the dependence on the angular law of scattering this formula 
may be found in the literature. Its most important special cases are 
obtained for t = const (Pseudo-Maxwellian distribution) and 
t = const/c (Druyvcsteyn distribution). 

The derivation of (65) should be completed by a proof that indeed 
h 2 (c) is small compared to fti(c). This is not true for all values of c; on 
the contrary, the argument below shows that near the origin where 
c ~ o,t(c), ho(c) is actually comparable to lh(c). For our purposes, however, 
it is sufficient if it is true in the overwhelming majority of cases. As the 
proof applies equally to all h,'s we will run it in this manner. Assuming 
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provisionally that 

h,+i <3C h„ <£ hy-i 

we see that we can drop the terms in h,+i in the system (47). The integral 
in (47) is evaluated in zero order with the help of (61) and (62); it 
becomes then 



^4-1 f P, (cos x)n(x) sin x rfx 

T(C) A Jo 

it 
The formula is 



r(«) 

This means that we may solve explicitly for h y (c) in terms of k-i(c). 



dhy-iic) v — 1 , , x 
3 1 n—i\fi) 

h (c) » — - - - (66) 

KM ~ 2v _ x / i - p y (cos X ) \ 

\ ar(c) I 

To estimate the order of magnitude of this we may neglect P, (cos x) as 
compared to 1 and assume v large. The operation in the numerator will 
lead to two kinds of terms: some of the form 

- h-M 
c 



and others of the type 

m c 



-- k-i(c) 



M {{ar{c)Y- 

coming from differentiation of an exponent of the type (65). Now we 
find from (65) that the overwhelming majority of particles have speeds 
c which in order of magnitude satisfy 

c - (M/m) 1/2 ar(c) (67) 

because this is the range within which the exponent remains comparable 
to 1. Applying this to the two types of terms arising from (66) we find 
for them in order of magnitude 



var(c) - K-i(c) r*t v I —J K-i(c) 



and 

,1/2 



^FiWb^ w ~© m ^ w 
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If we substitute this into (66) we see that the h p 's decrease as (m/M )* ", 
with a possible v\ slowing up the final convergence. In any case h 2 (c) 
comes out small compared to hi(c) which is all that is needed to make 
equation (65) approximately correct. 

While the case of small m/M is generally known it appears to be other- 
wise for large m/M. The intuitive basis for the solution of this case is 
the fact observable from (52), (57) and (59) that (c) increases indefinitely 
with m/M, but that the relative deviation from the mean decreases so 
that the distribution function approaches a 5-function. The structure of 
this limiting function may be explored, starting directly from (40), 
because in the limit of large m/M the sphere of integration shrinks as 
may be verified from (37). This makes it possible to replace the integral 
in (40) by differential terms. This becomes clearer if (40) is written in 
the form 

dh(c) 1 f* x . . 1 f" , //c'Y He') h(c)\ ,„. 

Let us call the inner average $. It exhibits the differential properties 
discussed earlier: the curly bracket is the difference of two terms which 
are almost identical. Hence we approximate the value of ^by expanding 
the slowly varying terms to first order in c' — c, while the fast varying 
h(c r ) will be expanded to square terms. This expansion is obviously 
permissible for everything except the rapidly varying function /i(c'). 
For h(c') itself no justification can be offered except success. By proceed- 
ing to square terms in this expansion we mitigate any possible error com- 
mitted, but it is quite possible that structural details are lost in the 
procedure. 

The development of $ is straightforward. We proceed as follows 

< " I f *• WSJ £) - ?fe) + (0 TO ^ ~ ^A 

~ m {©' W) - A + ^J'K f "" IMc,) - " (c)i 

The expansion of the first term involves only (41) which to this order 
reads 

d — cfHc — (1 — cos x) 
m 

This formula does not contain the azimuth w which therefore disappears 
trivially. In the second term on the contrary we are dealing with a 
vectorial difference involving all three polar coordinates c', k, oj of c\ 
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If we set 



C{ = c' sin k cos w 



c n = c sin k sin co 



Cf = C COS K 



we find 

ft(C) - fc(0 = — c £ + — c, + — ( Cf c) + 2 ^ « 

, i d 2 /t(c) , 2 i e*h(c) ,, _ )2 d^m , , 

, 3 2 /i(c) ,, / . , a 2 /i(c) /, / s , t 
With the formulas given, integration over a is elementary. We find 

, M , 1(c) i(cr i _ i \ + i (^ <„< co S . - c) 

\\c/ t(c') t(c)J t(c) 1 ac r 
^ 4r(c) \ d^ w ) 2r (c) dee 

All coefficients are to be evaluated only to the lowest non-vanishing order 
in M/m. From the equations (41) and (42) we get 

c' cos k — c « (M/m) c(l - cos x) 
c' 2 sin 2 k W (M/mfc sin 2 x 

The first term in J" is simplified further by introduction of the parameter 
a used in the dimensional analysis of Section ID. According to that 
section we have that 

dtnr{c) m _ 1 + a ( 69 ) 

d tnc 

We can generalize the original definition for any r(c) by the above 
equation, where a is now a function of c. Eliminating also the tempo- 
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rary device of a £, ij, f coordinate system we find for t <(: 

a M (a < \\/i \ ^( c ) . M (! _ cos x) ( dh , dh 

m T {c) m t(c) { dc x ^ dcv 

, ti\ , (M\ 2 sin^x fafa £h M[ 
' Be] \m) 4t(c) \dc 2 + dc u 2 + des 2 ) 

, /^Y 2(1 - cos x) 2 - sin 2 x / a 6 d ' 

W 4t(c) " ° r ^ + * dc~ y + Ca a^ 

c x dh Cy dh c z dh 
c dc x c dCy c dc z 

If the last term is evaluated exactly terms of the order (M/m) 2 are 
added to the first derivative terms in h(c). They are obviously negligible. 
Finally integration over x yields the following form for equation (68) 

dk(c) 



c) , ,> M l\ — cos x\ r . / ,,,/ » 



dc z m \ t(c) 

M / l - cos x \ ^ d* 1 /MV / sin 2 x\ 2 v a 2 /i , . 

, 1 (M\ / (l - cos x)(l - 3 cos x)\ ^ d 2 h 

~ t l\mJ\ r(c) /«**55K 

When equation (70) is considered up to linear terms in M/m it yields 
a 5-function about the drift velocity (c,) which results from the implicit 
equation 

The 5-function takes here the aspect of a non integrable function which 
in a special case can be seen to equal 

*« {c; + C J + (* - <c» 2 r 3/2 

When normalization is imposed on such a function it is made to vanish 
everywhere except at the point corresponding to the drift velocity. 

The square terms in M/m are necessary to gain information about the 
functional form of h(c). Since the region in velocity space in which h is 
appreciable is still small we may take a as constant in that region. A 
further simplication results from order of magnitude considerations 
on c: 

c 2 ~ cl ~ (c,f » cl ~ d ~ (cz - <c,» 2 
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In working out the details we see that division by 

M /l — cos x\ 
m \ t(c) / 

puts the first two terms on the right hand side of (70) in a simple form. 
The coefficient of the field term becomes then 

am am ({c^ 



\ t(c) / \ t{c) /«-<«.) 

The first factor is just (c«> by the identity (71). For the second factor we 
have up to linear small terms 

c - Vcl + c 2 v + cl = Vcl + cl + ((c) + c z - < C ,» 2 



and hence for the coefficient of the field term 

(c z ) (feJY " « (c z ) - (1 - a)(c z - (c z )) 

After division by 

il//l — cos x\ 
m\ t(c) /, 

the square terms still contain another small factor M/m; it appears 
sufficient, therefore, to keep only the leading terms which are the ones 
containing (c,) 2 as factor. All these terms are multiplied with the ratio 
of two angular averages over dx; these may be taken as independent of 
c to a good approximation. Equation (70) thus takes the form 

(4 - a)h(c) + c x — — + Cy — — + (2 - a)(c z - <c z » — — 
ac x oCy vvt 

/ sin 2 x \ 
\M \ t / / N 2 jd 2 h(c) d 2 h(c) \ 
h 4 m /I - cos x\ W I ^ + 64 / (72) 

— / 

/ (l - COS X)\ 

i 1 M \ T 1 1 V-^l = o 

+ 2m/ (l- cosxV ^ del 

\ ; / 
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The equation can be solved explicitly in Cartesian coordinates by the 
method of separation of variables. The result is 



h(c) = exp 



/ l — cos x \ 

—2 — ^ T ^ c * ~*~ c " 
M /sm2_x\ <c 2 ) 2 

\ T / 

/ l — cos x \ 

- (2 - a ) — -^ T ' (c ' "" ^ 



(73) 



M / (l - cos x ) 2 \ <c,) 2 
\ r I 

This is a Maxwellian distribution with elliptic distortion and shifted 
origin, that is, the type shown in Fig. 1 (b). 

The result (73) indicates the main features of the solution for heavy 
ions. Because of the neglect of derivatives higher than the second in 
/i(c') it is not certain that (73) is correct in all details, even in the limit 
of very large m/M. 

IID. THE CASE OF EQUAL MASSES; IONS TRAVELLING IN THE PARENT GAS 

The developments of the previous section show that if the ion mass is 
either large or small in comparison to the molecular mass, analytical 
methods can be applied successfully to determine the velocity distribu- 
tion of the ions. No such possibility' was found for the mass ratio unity, 
which one would judge to be of particular interest because it applies to 
ions travelling in the parent gas. There exist isolated fragments of such 
an analytical theory; for instance, if we assume isotropic scattering in 
(46), that is n (x) = 1, then the zeroth equation (46) becomes explicitly 
integrable and yields 

He) = - | ar(c) *$ ■ (74) 

This is a curious reversal of the differential relationship (63) derived for 
electrons and implies a rather strong condition on the structure of 
hi(c). However, I have not been able to consolidate these fragments into 
something which can be used successfully in computation. The high 
field distribution function for mass ratio unity appears, however, suffi- 
ciently interesting to warrant the use of other methods. 

A numerical determination of the velocity distribution function was 
undertaken in cooperation with R. W. Hamming by the so-called 
Monte Carlo method. The Monte Carlo method is a way of gaining 
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statistical information about a system by following an individual member 
through a large number of random processes. The result of such a pro- 
cedure is knowledge about one member of the assembly for a long period 
of time. Time averages of various kinds can be obtained from such data; 
these time averages are then set equal to instantaneous averages over the 
assembly, in accordance with ergodic theory. In our case, an ion was 
followed through 10,000 collisions. On an average, the collisions were 
isotropic in the center of mass system (n(x) = 1) and obeyed a mean 
free time condition r = const. Actually, both the free time and the scat- 
tering angles varied from collision to collision; the angles varied in a 
random fashion over a unit sphere and t was random within an ex- 
ponential distribution. 

A Monte Carlo calculation of this type consists of three parts. In the 
first part the random numbers having the required distributions are 
obtained and recorded. In the present problem there were three such 
random numbers required for each collision, namely a time and two 
angles. These numbers were placed on 10,000 IBM cards, along with 
suitable identification. In the second part a calculating machine simulates 
the successive collisions and keeps a record of the initial and final veloci- 
ties for each one. The third part consists in analyzing statistically the 
numerical material accumulated in the second. For the first part of the 
calculation particular values must be chosen for the acceleration a and 
the mean free time r. These values were 

a = 1 

r = log w e = 0.43429 

However, the dimensional analysis of Section ID shows us at this point 
that these two constants enter into the problem only through their 
product ar which scales all velocities. It is therefore convenient at the 
statistical stage to remove these factors and to analyze the results in 
terms of a dimensionless variable which by (26c) we take in the form 

w = ■£ (75) 

ar 

In view of the a priori information for mean free time problems which 
is gathered in Section IIB we can use the statistical data from the 
Monte Carlo calculation in two ways. "We may (a) check the numerical 
computation itself or (b) gain new information not available otherwise. 

(a) The check of the numerical calculation by statistical analysis 
proceeds as follows. From deductive reasoning we have obtained the 
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averages (52), (54) and (57) for c z , c\ and c\ . These formulas ought 
to be verified in the Monte Carlo calculation. This is indeed approxi- 
mately true. A sampling covering 9492 out of the 10,000 collisions gives 

M 
W) 

<-;> 

The agreement is essentially there but deviations are noticeable. In 
judging these we have to realize that fluctuations are quite large in this 
problem. For instance if the calculation is broken down into ten runs of 
approximately 1,000 events each one finds the following time averages 
for the partial runs: 

<«-«> («£) <«5) 



by Monte Carlo 


by deduction 


1.912 


2 


0.801 


8 

- =0.889 

>' 


5.165 


56 

— = 6.222 



1.821 


0.837 


4.453 


1.975 


0.748 


5.531 


1.915 


0.785 


5.132 


1.954 


0.829 


5.441 


1.868 


0.731 


4.794 


2.003 


0.807 


5.581 


1.766 


0.793 


4.811 


1.962 


0.827 


5.360 


2.003 


0.901 


5.471 


1.914 


0.727 


5.200 



predicted 2.000 0.889 6.222 

Among these runs there are some having averages higher than the pre- 
dicted values, but the data clearly show that the Monte Carlo averages 
are generally lower. In the search for reasons it was first felt that perhaps 
the desired mean value for r is not actually reached, perhaps through 
systematic errors introduced by the operator when rejecting certain 
runs. This seems indeed to be the case. The mean free time obtained 
from the 9492 runs mentioned above is 

t = 0.4269 

•which is slightly low. Indeed it is observed that the runs with high r 
were particularly troublesome in the calculation and were preferably 
rejected by the operator. It seems doubtful however that this error 
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could account for the entire discrepancy, particularly in the mean 
squares, although it must be emphasized that the runs with high r make 
a more than proportional contribution to the total average. The angles 
of scattering have not been subjected to a similar analysis so that we 
cannot make a statement whether the aimed at isotropy in the law of 
scattering was realized or not. We conclude therefore by saying that 
while the Monte Carlo calculation gives results in general agreement 
with the deductive theory there are small but noticeable systematic 
errors in it whose origin is only partly explained. Similar errors must 
exist in the new results which cannot be compared with theoretical 
predictions. 

(b) In this part we will discuss the velocity distribution function which 
may be constructed from the Monte Carlo results. In constructing such 
a function we make use of the fact that, between collisions, the velocity 
is accelerated at a uniform rate. Thus, in each period between two 
collisions, the velocity vector traces out a straight line parallel to the w 2 
axis covering equal distances in equal times. The Monte Carlo calculation 
furnishes us with a number of such straight lines as shown in Fig. 1 1 . The 
density of these straight line tracks in velocity space is the velocity 
distribution function. The actual procedure used to obtain it was to lay 
a grid with a mesh of 0.23 in a half-plane with coordinates w t and 
w p = -\/wl + w\ and to count the number of lines crossing each hori- 
zontal square edge. When the resultant count is converted to density 
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% = VW X 2 + Wj 



Fig. 11 — Straight line pattern in the w t - w, half plane from which the veloc- 
ity distribution is constructed; the Monte Carlo calculation furnishes the initial 
and final velocities (dots and rings). 



MOTION OP GASEOUS IONS IN STRONG ELECTRIC FIELDS 



215 



and normalized to 1 we get a distribution in the w z — w p half plane which 
is shown in Fig. 12. Division by 2mv t will transform it into a conventional 
distribution function in velocity space; a plot of this function is shown 
in Fig. 13. What distinguishes this distribution function from functions 
previously proposed is the elongated probability contours. This feature 
is not unexpected in view of the unequal energy partition apparent in 
the equations (58) and (60). 

The probability contours shown in Figs. 12 and 13 give a reliable general 
picture but we must not expect from them fine detail. Indeed we will 
now prove that the distribution function is infinite along the entire 
positive Cj-axis, a feature which is not obvious from inspection of Fig. 13. 
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Fig. 12 — Motion of ions through the parent gas in a high field; distribution of 
velocities in the w, — w, half plane resulting from the Monte Carlo calculation. 
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W x 

Fig. 13 — Velocity distribution function of ions moving through the parent 
gas in a high field; contours constructed from the Monte Carlo results of Fig. 12. 

A simple physiccal proof of this statement goes as follows. Suppose an 
ion and a molecule make a collision which is almost central, but has a 
small impact parameter b. The collision will bring the ion almost to rest 
because the atom was originally at rest by hypothesis. Because the col- 
lision was not quite central, however, the ion will have a small residual 
velocity c f at right angles to its original velocity c, . For any reasonable 
law of scattering this quantity c f will be proportional to c< and to b. 

C/ cc Cfb 

The probability for a value of b between b and b + db is proportional 
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to b db. Thus even if all d's were equally probable the probability for 
c/ would vary as c f dc f . Actually very small Ci's may be specially probable 
as the theorem states and this fact may or may not increase the prob- 
ability for small cj. This means that P(c/)dc/ probably varies as c/dc/, and 
may perhaps even contain a smaller power of c f . When such a probability 
function is plotted in velocity space it will vary as \/c f . Thus we know 
that the distribution function <p(c) for ions immediately following a col- 
lision has a singularity at the origin at least as 1/c/ . The actual distri- 
bution function h(c) is derived from this one by spreading each point 
out in the forward direction as shown in Fig. 1 1 . For the mean free time 
case we can write this out explicitly in the form 

1 r w * 

/i(c) = - / <p(c - at) e r dt (76) 

T JO 

For the case of a mean free path or other laws the formula is more awk- 

i 
ward but they all differ from the above only in replacing e f by a more 

complicated weight function. The singularity of h arises out of the singu- 
larity of <p which contains at least a factor l/\/& + Cy + (c 3 — at) 2 . 
Along the c,-axis, this become a factor l/(c z — at); this factor makes the 
integral diverge for all positive c z ; as we approach the origin from nega- 
tive CgS the distribution function will become infinite at least as in 1/c, . 

The reasoning given is intrinsicly classical because of the use of "in- 
finitely small" impact parameters. We should not hasten to conclude, 
however, that the quantization of the angular momentum will necessarily 
remove the singularity. Indeed we know that the only mechanical 
information which has to be put into the Boltzmann equation (31) is the 
differential cross-section for scattering. If this quantity does not differ 
essentially in the 180° direction from a classical cross section then it 
will not modify the conclusion we have reached. 

Conclusions which are more informative, but less "anschaulich" may 
be obtained from a study of Boltzmann's equation either in its closed 
form (34) or (40) or its "Legendre" form (46) or (47). In view of the 
proof given we will give only an outline of the reasoning. First we can 
remove the s€;cond term in (40) by the substitution 

dc z 
ot(c). 

The exponential is easily seen to be always positive and finite for finite c. 
The Boltzmann equation then takes the form 



h(c) = exp[-£' "' |/ c 



f f c ' dc z 1 dh* 1 , . , n 

ex P ~ / rrzs \ a ir = - ' ( an lnte g r ai) 

L Jo ar(c)J dc z c 
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For equal masses the integral on the right runs over a plane in velocity 
space. Its integrand is always positive; hence the integral can never 
vanish and is always positive; it is conceivable that it could be infinite 
for a special position provided the infinity is integrable. Such an infinity 
would only make matters worse. At any rate the equation shows that 
h*(c) increases monotonely with increasing c, . When we approach the 
origin from negative c z we get a logarithmic divergence (or worse). As 
the function h*(c) can nowhere decrease with increasing c* the infinity 
along the positive c, axis is confirmed and its logarithmic nature is made 

very likely. 

Information obtainable from equation (46) confirms this conclusion. 
Lowest powers in the entire recursion system can be made to cancel by 
assuming that for small c 

ho ~ — Alnw 

hi ~ Bi i > 

with suitable relationships existing between these quantities. 

A defect of all three approaches is that they give no information 
concerning the nature of the infinity for c z > 0. One is tempted to 
conclude from Fig. 13 that it cannot be very strong. Something like a 
singularity is discernible at the origin, particularly if the contour 0.1 is 
drawn back to cut the w a -axis at a negative value; this is perfectly com- 
patible with the available information. For large positive w, , on the other 
hand, the picture almost contradicts the theorem just proved. One con- 
cludes from this that the singularity, for large c, , becomes a weak and 
narrow ridge rising more or less abruptly in an otherwise well behaved 
function. 

HE. THE CASE OF EQUAL MASSES; A NEW COMPUTATIONAL PROCEDURE 

The foregoing sections have accumulated substantial evidence that 
there are many analytical details involved when one discusses the 
structure of a velocity distribution function. These details are of little 
interest to the experimenter who may want nothing but a formula for 
the drift velocity or the average energy. In view of this situation it 
appears very desirable to find a method whereby such quantities can be 
derived directly and accurately from the Boltzmann equation without a 
full knowledge of the entire distribution. 

Maxwell's original work shows us how to achieve this for molecules 
obeying the mean free time condition of Section IIB. In the following, 
a general method is described which will permit determination of such 
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averages for an arbitrary law of force between the ions and the gas 
molecules, and an arbitrary mass ratio. The application will be limited 
to the case of the mass ratio 1 whose study was begun in the preceding 
section. 

The basis of the method is an observation on the equation system (46) 
or (47), which is the form taken by the Boltzmann equation after in- 
serting the Legendre decomposition (43). It would appear at first sight 
that these recursion relations are of such a structure that an arbitrary 
function ho(c) could be substituted into the "zeroth" equation and that 
the relations would then successively determine hi , hi , h 3 • • • . Upon 
closer inpsection this is found not to be the case. Suppose we have 
obtained somehow functions h , hi , ha • • • h n and we are trying to use 
the nth equation to determine h n+ i . This equation is of the form 

— ? h h n+ i = known material (77) 

dc c 

We solve for h n+ \ by multiplying with c n+2 and integrating. This gives 

c n+2 h n+ i(c) = I (known material) dc 

The left-hand side is of such a structure that it must vanish both for 
c = and c = <x> . It follows that the right-hand integral when taken 
between the limits and °o must equal zero. This condition is indeed 
obeyed for any ho(c) when n = 0. The integral condition reads in this 
case 

- / c dc n(x) sin x dx — r-rt ( - - / — ^f c l dc = 
2 Jo Jo t(c ) \c / Jo t{c) 

If we invert the order of integration in the double integral, then intro- 
duce c' as variable of integration by equation (41) and finally invert 
again this becomes 

This equation is trivally obeyed because the square bracket vanishes in 
virtue of the definition of IT(x). For values of n higher than 0, the in- 
tegrability condition deduced from (77) is not generally obeyed for any 
function ho(c). Such a statement may be proved by examples; these 
examples will arise in the course of the calculations to follow. Thus 
we find that except in the passage from ho(c) to hi(c), the recursion system 
is such that at each stage it imposes a condition upon the h,'s already 
determined if the new h n+ i(c) is to exist at all. With such an infinity of 



220 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1953 

conditions one can improve indefinitely an initial trial function assumed 

for ho(c). 

The integrability conditions whose general structure is thus indicated 
have actually already been written down. They are the equations (49) 
for the special case s = v. Generally speaking, the relations (49) are also 
of the recursion type, permitting us to start with arbitrary averages 
(c), and computing successively (c"Pi (cos t?)) etc. At every stage, how- 
ever, there is the exception mentioned: the equation for which s = v 
has no third member, and therefore it imposes a condition upon averages 
already known from the previous equations. We shall refer to this type 
of equation as a "truncated" relation. 

It is reasonable to assume that 1/t(c) can be developed into a power 
series in c because it equals the known constant polarization value for 
c = 0. If this can be assumed then each truncated relation s = v is 
equivalent to a unique relation among velocity averages involving 
ho(c) only. One obtains this relation by applying to each member in the 
truncated relation its own recursion formula and repeating this process 
until v is brought down to zero. This process will never lead into another 
truncated relation s' = v' because at each step s' increases by at least 
two units with respect to v' '. 

In order to test the method for a known case, it will be applied first 
to the case of constant mean free time. This case is adequately described 
by the theoretical treatment of Section IIB and the Monte Carlo 
calculation of Section IID. We have seen that the equations (49) reduce 
in this case to the form (51) which dovetails as shown in Fig. 9; this 
dovetailing leads to explicit values for certain averages as shown in 
Fig. 10. A "computational method" is only needed when one tries to get 
an average outside this selected list. In the present case the reduction of 
the truncated relations to a condition on h(w) is particularly simple as 
is seen from Fig. 9. A singular relation which starts out as between 
(c~ l P r -i (cos #)) and (cP, (cos t?)) actually yields the numerical value 
of the latter because the former has been obtained numerically in a previ- 
ous stage. This numerical value yields in combination with previous 
information (c +1 P,-i (cos #)), (c" +2 P F _ 2 (cos &)) etc and finally (c"). 
Thus we end up with the set of even moments of ho(c) which may be 
used in succession to determine ho(c) more and more closely. There is no 
guarantee that this procedure converges mathematically, since the 
general theorems usually require the knowledge of all integer moments. 

20 Shohat, J. A., and J. D. Tamarkin, The Problem of Moments. Am. Math. 
Soc, 1943. The original three-dimensional formulation appears a little more 
favorable for a proof because, in this case, we know indeed all integer moments. 
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The justification for the method rests therefore on an empirical basis at 
this point. 

Assuming isotropic scattering, as in the "Monte Carlo" calculation 
we express our results in terms of the dimensionless variable w defined 
in (75). The equation system (51) becomes then 

(2, - 1)(1 - </..,(x)»<«, r) - 

(78) 
«,{> + «+ l)(s - 1,-v - 1) f (» + 1)(„ - s )(s - 1, v + 1) 

where the abbreviation (s, v) has been introduced for (w*P, (cos #)) 
and the quantities (/«,,(x)) are simple numbers computable from (48b) 
and the assumption of isotropic scattering. The first truncated relation 
is s = v — 1. It yields 

(1, 1) = 2 

Reducing it with the relation (78) for which s = 2, v = we get 

(2, 0) = 8 (79) 

The next truncated relation is s = v = 2, which yields 





**-¥ 


and the reduction gives 






<3, 1> = | 




(4, 0) = 184 


Similarly in the next stage 






(3, 3) = 1H§ 




1A OX 18112 

<4 ' 2> " 133 




,_ tX 421600 
(0 ' 1>_ 399 




/r „ x 3372800 
< 6 ' °> 5J5T 



(80) 



(81) 

As an example of an average which cannot be had explicitly we may 
take the mean absolute value of the speed, that is (1, 0). We find this 
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value by picking a sequence of trial functions for ho(w) with the ap- 
propriate number of parameters and imposing successively (79), (80) and 
(81) upon this sequence; this leads us to a sequence of values for (w) 
which can then be examined. In such a procedure careful consideration 
of the trial functions is an important element. The following information 
is available. It was proved in Section IID that hv(w) is logarithmically 
infinite at the origin. At infinity, on the other hand, h(w) falls as e " 
times some power of w. One way to check this is to drop the terms con- 
taining 1/c as factor in (46) ; the solution of the recursion system becomes 
then 

h,(w) ~ (2v + l)e' w w k 

where k is some unknown exponent. Armed with this fore-knowledge, 
we shall use the following sequence of trial function for ho(w) 

ho{w) = pEi(w) + qK (w) + re~ w + swKi(w) (82) 

where 

.00 — u 

Ei(w) = I — du 

Jw u 

and K (w), Ki(w) are the modified Hankel functions of order zero 
and 1. 
•We find in zeroth approximation from normalization only 



(83a) 



v m = \ 


q m = r m 




(w) m = 2.2500 


in first approximation, using (79) 




p m - 8 




*'" = t 




(w) {i) = 2.3818 



r (1) = S (1) - 

(83b) 



J1 This definition, which is in accord with the tables of Jahnke-Emde, differs 
from the usual one by a factor 2/sr. This change is suggested by Watson, Bessel 
Functions, p. 79, and proves convenient in the following. 
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in second approximation, using (79) and (80) 



p (2) = 1 

P 12 




(2) 32 
1 45 


8 l% 


r (2 > - _ 1 
20 




(w)™ = 2.3858 




and in third approximation, using' (79), and (80) and (81) 


o) _ 3079 
P 7980 




(3 ) _ 202544 
q 209475 




( 3) 2507 
18620 




( t) 3152 
209475 




(w) 3 = 2.3864 





(2) = 



(83c) 



(83d) 

Appearances indicate strongly that the sequence (83) for (w) approaches 
a limit which one would guess to be 

(w) = 2.3865 (84) 

More evidence that the conclusion drawn is correct can be obtained 
by using the set of trial functions 

ho(w) = pK (w) + qe~ w + rwe~ w 

We find then the following sequence of values for (w). 

(w) m = 2.546 (w) {1) = 2.395 <w) (2) = 2.388 

This descending sequence confirms (84) by approaching this same value 
from above. 

Further evidence for the correctness of the procedure can be obtained 
by deriving a function ho(w) from the Monte Carlo function h(w) dis- 
cussed in Section IID and comparing it with our trial function. The 
function was constructed by covering Fig. 13 with a grid of concentric 
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ho(w) = 2x / h(w) d (cos &) 

J 6=0 



circles and horizontal lines and replacing the integration 

..>=■»■ 

by a summation over grid points. The function ht(w) so obtained is 
compared in Table I with h^iw), ho ll \w) and k> l2) (w) as denned by 
(82) and the numbers following. We observe that the first approximation 

Table I 

Comparison of the Monte Carlo ht(w) for h (w) with successive approximations 
obtained by the new method. 



w 


*?(«>) 


WHw) 


*$*>(») 


itfHw) 







00 


co 


OO 


0.5 


0.74 


0.8397 


0.7281 


0.7144 


1 


0.29 


0.3291 


0.3019 


0.3002 


1.5 


0.15 


0.1500 


0.1438 


0.1440 


2 


0.081 


0.0734 


0.0730 


0.0733 


2.5 


0.0412 


0.0374 


0.0384 


0.0387 


3 


0.0199 


0.0196 


0.0207 


0.0208 


3.5 


0.0118 


0.0105 


0.0114 


0.0114 


4 


0.0063 


0.0057 


0.0063 


0.0063 


4.5 


0.0030 


0.0031 


0.0035 


0.0036 


5 


0.0014 


0.0017 


0.0020 


0.0020 


6 


0.0004 


0.0005 


0.0007 


0.0006 


7 


0.0001 


0.0002 


0.0002 


0.0002 



is an improvement over the zeroth one, while the second one makes little 
difference, considering the accuracy to which fc o («0 is given. In indi- 
vidual cases the sequence drifts away from Iit(w) ; this is not surprising 
because the latter function is very rough; this is to be expected from its 
mode of derivation. 

The application of this method to the hard sphere model of ion-atom 
collisions offers no new feature of principle. The actual working out of 
results is somewhat more complicated, mainly because the connection 
diagram for the recursion system (49) is more involved. According to 
equation (26b) the dimensionless variable to be used in this work is 

(85) 



\/aX 
We denote its averages (w'P, (cos &)) by (s, v) as previously. The equa- 
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tion system (49) then takes the form 

( 2 „ + 1)(1 - (I.„( X )»< 8 + 1, v) = 

= v(v + 8 + l)<s - l r - 1) + ( r + l)( s - y )(s -!,„+!} 
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(86) 



The numbers (I,,v(x)) were already discussed in connection with the 
system (78). What distinguishes (86) from (78) is the way in which the 
variables are connected; the new connection diagram which replaces Fig. 9 
is shown in Fig. 14. The truncated relations no longer dovetail into each 
other as they did before. Only the first stage proceeds in a similar way, 
yielding explicit expressions for (2, 1) and (4, 0). In the next stage we 




Fig. 14 — Interconnection established by the Boltzmann equation among the 
averages (c'P, (cos #)); case of constant mean free path. 
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start out with a relation between (1, 1) and (3, 2). By the use of regular 
recursion formulas we can successively transform this into a relation 
between (3, 0) and (3, 2), then (3, 0) and (5, 1) and finally between 
(3, 0) and (7, 0). Here we have for the first time the normal situation in 
which we do not get the actual value of a moment of /io(c) but only a 
relation between two or more of such moments; the reason for this is 
that the system fails to connect up with (0, 0) which equals unity a 
priori. A similar situation prevails for the next truncated relation; it is 
originally a relation between (2, 2 > and (4, 3 ) and is finally reduced to one 
between (2, 0), (6, 0) and (10, 0). Similarly, the next truncated relation 
reduces to a relation between (5, 0), (9, 0) and (13, 0) and so forth. The 
first three of these reduced relations come out to be 

(w*) = 10 (87) 

3<u> 7 > = 112<w 3 > (88) 

** (w<) = 27^) + il (w 10 ) (89) 

These formulas will now be imposed upon a sequence of trial functions 
for h (w) suitably chosen. Again, we may make use of the information 
of Section IID, according to which ho(w) is logarithmically singular at 
the origin. For large w we proceed as previously from (46) leaving off the 
terms of 1/c. We get then 

K{w) ~ (2v 4- l)e~* "V 
This suggests the following trial function for h{w) 

ho{w) = pEi{W) + qKoiW) + w H *' + 8W*Ki(hu?) (90) 
The best zero order approximation is actually obtained by the function 
K (W)- We find 

«'" = 2T^) = L04605 p(0> = r(W = s<0> = ° 



r = s 



In first order we get, using (87) 

p (1) = -0.46543 

q {i) = 1.45285 
In the second order, using (87) and (88) 

p (2) = -0.80856 

? (2) = 1.88127 s <2) = 
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In third order, using (87) and (88) and (89) 

p (3) = -1.15071 

q (3) = 2.37034 

s (3) = 0.02062 

r m = -0.29016 

These successive approximations lead to the following sequence for 
the drift velocity (w cos t?) 

(w cos #) (0) = 1.04605 (91a) 

<wcos#) (1) = 1.14256 (91b) 

(w cos #> (2) = 1.14616 (91c) 

(w cos tf> (3) = 1.14661 (91d) 

We conclude from this sequence that 

(w cos #) = 1.1467 (92) 

In addition to the drift velocity there is some interest in the energy 
and the energy partition. For the energy the following numbers are 
obtained 

<w 2 ) (0) = 2.1884 (93a) 

<w 2 ) (1) = 2.3395 (93b) 

(w 2 ) l2) = 2.3511 (93c) 

<w 2 > <3) = 2.3531 (93d) 
giving 

(w*> = 2.353 (94) 

A zero order value for (w 2 cos 2 t?) cannot be said to exist because the 
first truncated relation is the condition that a distribution function 
Jh(w) exists at all. Thus, we can get only three numbers in a sequence 
approximating (w 2 cos 2 #) 

(w 2 cos 2 #) (1) = 1.8005 (95a) 

<w 2 cos 2 t>> (2) = 1.7696 (95b) 

(w 2 cos 2 #> (3) = 1.7685 (95c) 



' 
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giving 

(w 2 cos 2 &) = 1.768 (96) 

We can understand the results (92), (94) and (96) by giving the frac- 
tion of the total energy in ordered motion and the fraction of the energy 
in motion along the z-direction. We find for the first ratio 



<»"»»> = 0.559 (97) 



and for the second 



(w cos fl) _ 0751 (98) 

The ratio (97) equals 0.5000 for all mean free time models; the ratio 
(98) is 0.778 for the mean free time model with isotropic scattering. 
Thus, the deviations from the earlier results are not drastic. However, 
in certain derived relations the difference is more noticeable. For instance, 
a good measure of the anistropy of the diffusion process is furnished by 
the ratio of the random energy along the field to the energy at right 
angles. 22 From (97) and (98) we find for this number 

(W 2 COS 2 g) - (W COS ff) 2 m 1 ^ (gg) 

KM - (w^cos 2 ^)) 

For the mean free time case this number equals 2.50. Hershey in his 
work assumes this number to be 1.000. 

A comprehensive list of velocity averages is attached in Table II. As 
a comment I may add that the obvious mode of constructing such a 
table, namely by computing the column v = from (90) and then using 
the recursion system (86) for the others, runs into some difficulty. First 
of all, a series of cancellations reduces the accuracy as v increases; 
finally, at the positions marked "impossible" we find the missing third 
members of the truncated relations. These elements cannot be com- 
puted by recursion at all, but would require an explicit solution of the 
equation system (47) for h,+i(c). In the table, this more arduous path is 
not followed. Instead, the recursion method is used for the numbers in 
italic type and a few numbers are added by extrapolation. The numbers 
so obtained will be needed in Section IVB. 

The calculations on the hard sphere model are immediately applicable 
to the experimental data of Figs. 3 to 7, which exhibit the drift velocity of 

22 See below, equations (147) and (165). 
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the noble gas ions in the parent gas as functions of the parameter a/N. 
These data have a high field range in which the drift velocity varies as 
the square root of a/N. This is the variation for a model with constant 
mean free path, as seen from the dimensional formula (26b). It was 
indicated furthermore in the Section IA that we have good reason to 
think of the scattering between an ion and an atom as nearly isotropic. 23 
These two features characterize uniquely the hard sphere model whose 
treatment we have just completed. To the extent that they are verified 

Table II 

Dimensionless high-field velocity averages (s, v) for the hard sphere model 
and mass ratio unity. 



V 





l 


2 


3 


4 


8 



1.0000 


0.7845 


impossible 






1 


1.S92S 


1.1467 


0.8022 


impossible 




2 


2.868b 


2.0000 


1.4769 


0.990 


impossible 


3 


4.5868 


3.9858 


3.0578 


2.134 




4 


10.0000 


8.8853 


6.992 


5.0602 


3.474 


5 


28.912 


21.405 


17.330 


12.84 




6 


61.847 


55.97 


46.177 


35.36 




7 


171.241 


156.3 


130.91 






8 


503.7 


462.81 








9 


1568 


1445 








10 


5090.9 


4750 









the model is applicable to the experimental data. The formula to apply 
is (92) in combination with (85): 



<ft>- 



-™j/£ 



(100) 



In the logarithmic plot of (c*) vs a/N the intercept of the straight line of 
slope 1/2 which fits the high field data thus equals 



log 



1.147 



The values for a which result from this are shown in Table III. For 
comparison are shown the corresponding atomic cross section as deter- 
mined from viscosity data. It is interesting to observe that the ratio of 

23 A quantitative discussion of this point for the polarization force will follow 
in Section IIIB. 

24 Landolt-Bornstein, 1950 edition, Vol. I, part 1, page 325. 
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the two retains very nearly the constant value 3 throughout the table. 
The fact that the ratio is substantially larger than unity is explained 
by the resonance feature of the ion-atom scattering process as discussed 
in Section IA. The fact that it is constant is perhaps an indication of 
the fact that both processes are governed by overlap conditions of 
essentially the same wave functions. 

I would like to point out in connection with the calculations of this 
section that the method developed is potentially of very wide applica- 
tion. One question that comes up, for instance, is whether a careful 
kinetics calculation is necessarily restricted to certain models or whether 
an ion-atom cross section known numerically could be used to derive 
therefrom kinetic properties. This is indeed possible. Suppose, for 
instance, that the cross section a(c) were available as a function of c for 
collision of He + -ions and He-atoms and suppose that this cross section 
were to satisfy the condition of isotropy n(x) = 1 to a good approxi- 

Table III 
Cross sections for ion-atom and atom-atom collisions for the noble gases. 



Gas 


ion-atom cross section X 10" cm' 


atom-atom cross section X lO" cm" 


He 


54 


15.0 


Ne 


65 


21.0 


A 


134 


42.0 


Kr 


157 


49 


Xe 


192 


67 



mation; we may then derive for this eventuality conditions on hn(c) 
which are more general, respectively, than (79) or (87), (80) or (88), (81) 
or (89). Since we are outrunning here the experimental evidence we shall 
limit ourselves to the derivation of the first of these relations. The first 
truncated relation is exactly (50a) which, for isotropy and equal masses, 
reads 

-SJ) = 2 (101) 

\(it(c)/ 

The reduction of this formula to a condition on k>(c) requires the rela- 
tionship v = of the set (46). This relation is always integrable to 
yield h{c) in terms of fto(c), as was pointed out early in this section. 
For the special circumstances assumed the integrated equation is 
equation (74) 

Hy) 



He) = 3 J 



ar(y) 



dy 



(102) 
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The elimination of hi(c) is achieved by forming the average (101) on the 
function (102). This leaves the required condition on ho(c); it may be 
given the following form 

The equations (79) and (87) are manifestly special cases of this more 
general relation. Adaptations of this procedure to other cases are clearly 
possible whenever the need arises. 

The calculations of this section are meant to suggest that it is possible 
to compute reliably average values from a Boltzmann equation without 
solving it completely. The method employed here for this purpose 
resembles a Ritz method in that it works with trial functions which must 
be guessed at, and like that method it is capable of indefinite improve- 
ment. The numerical results suggest strongly that we are converging 
toward a definite answer; however, a mathematical proof of this fact has 
not been presented. The method will be applied once more in the section 
on diffusion. 

Part III — Motion of Uniform Ion Streams in Intermediate Fields 

IIIA. A CONVOLUTION THEOREM 

Whenever we deal with the motion of a given type of charged particle 
in a gas of given composition, then there exists a wide range of densities 
n and N as discussed in Section IA in which the motion of these particles 
depends only on a/N and kT. For this range the motion is governed by 
equation (13). Since deriving that equation, all our efforts were dealing 
with the "high field" equation (34) or (40), in which the gas temperature 
is taken to be zero and the electric field often scales out, as in (26), (75) 
and (85). The accomplished solution of this restricted problem, together 
with the low field solutions available in the literature, brings us back to 
the more general equation (13) and the question what can be done with 
it. The topic of Part III so defined is definitely inferior in importance 
to the one in Part II. For we are studying here an intermediate range 
of variables which can be handled qualitatively, both in concept and 
practice, by some sort of interpolation between the high and low field 
regimes. For precise measurements, conditions can always be chosen so 
as to satisfy one or the other of the two extremes. For this reason the 
intermediate field case will only be pushed as far as it will go con- 
veniently, without appeal to numerical methods. 

In this Section IIIA we shall give a complete solution of the inter- 



(It 
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mediate field problem for the mean free time models discussed in Section 
IIB. This solution is achieved by the following theorem: Given the 
general equation (13) for constant mean free time 

d _M + /( c ) = _L [f M (CW)II(x) dSty dC (104) 

dC z 47T J J 

and the "high field" equation derived from it by setting the gas temperature 
equal zero 

a T *£@ + h(e) - J- [( «(C')A(cOn(x) dOy. dC ■ (105) 

and the Maxwellian equation derived from (104) by dropping the field 
term 

m(c) = i Jj M (C>(c')n(x) d$V dC (106) 

Z/ien </ie solution /(c) 0/ (10£) is i/ie convolution of the solution h(c) of 
(105) and the solution m(c) of (106): 

/(c) = j h(u)m(c - u) du (107) 

We carry through the proof by constructing explicitly the equation 
satisfied by the convolution. We replace the running variables c, c', C, 
C' in (105) by u, u', U, U' and multiply in m(c - u). We get 

a T *hp± m(c - ll) + h(u)m(c - u) = 
dUz 

= 4^ // ^ U 'W u ' )m(c " u > n ^u) dQ,. dJJ 

We now define /(c) by the relation (107), and integrate the above equation 
over u. The second member on the left comes out to be /(c). For the 
first member, we carry out an integration by parts: 

f a «5) m(c _ u)da= _/- Ku) am(o-u) du 

J du z J ou z 

= +f h (u) dMc - u)) 



= ± f h(u)m(c - u) 
dCt J 



dc 

_ a/(c) 

dc, 



dCz 

du 



du 
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For the right hand member we observe that we have the eightfold 
integration 

dn„- dU du, 

that is an integration over the collision angles and all final velocity 
components. By a general principle of kinetic theory 26 we can invert in 
this integration the final and the initial quantities and write 

d£l v > dU du = dQ, dU' du (108) 

This puts us in a position to eliminate the 5-function by integration. 
We find 



(IT 



^r + /(c) = h II ft(tt ' )w(c - u)n(x «) d ^ dW ( 109 ) 



with the side condition that u, U, u', U' form a quadruple of vectors 
in the sense discussed in Section IB for which in addition 

U' = 

If we substitute (107) into (104), denoting the dummy variable by u' 
instead of u, then the two equations (104) and (109) take on a very 
similar look. A proof of their identity hinges upon proving the identity 
of the integral terms : 

J h(u') du' f ro(c - u)n(xu) dO, 

(110) 
= J h(u) du JJ M(C')m(c - u')n(xc) da y dC 

The form of this relation suggests the assumption that the expressions 
are identical before integration over u'; this assumption is proved by 
the events below. The complicated function h(u) thus disappears from 
the problem. The other such function, namely n(x) disappears then also; 
for it is by assumption arbitrary, hence could be replaced by a 5-function 
for a fixed, but arbitrary x- The two sides of (110) must therefore be 
equal before we integrate over Xu or x c , and the two x's are to be taken 
equal and fixed. Defining angles as shown in the spherical diagram Fig. 
15 we thus get (110) in the form 

J to(c - tt) de = JJ M(C')m(c' - u') d<j> dC (Ilia) 



25 See Reference 4, Section 3.52. 
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This is to be true with the side conditions 

c = fixed (Hlb) 

u' = fixed (Ulc) 

U' = (Hid) 

Xc = Xu = X = fixed (Hie) 

Equation (111) is an identity involving only elementary functions. 
Thus the relation itself is in a sense elementary. Those who wish to 
believe it, may consider the theorem proved; for completeness, however, 
the proof of (111) will now follow. 

Call the left side of (Ilia) X, the right side Y. To determine X, we 
substitute from (7) and (Hid) 



u = 



m 



M -f- m 



u' + 



M 



M + m 



V 



with 



2 n.' 2 

■q = U 




Fig. 15 — Definition of the angles occurring in the proof of the convolution 
theorem. 
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because of (9). This yields with the angles as shown on Fig. 15 

X = / m(c — u) de = I — J exp — /3mc 2 

. ,2 M 2 + ra 2 + 2Mm cos x . 00 / m + M cos x .1 

- /Swim 71 r r - h 2Pmcu — — — cos ^ 

(M + m) 2 M + m J 

• / exp 2/3 irr ^- cu sin x sin d> cos e de 

Jo \_ M + m J 

The integral is evaluated by a formula known from the theory of Bessel 
functions 

[' e !C0Bt de = 2irh(z) (112) 

Jo 



and yields 

v 2 fa \3/2 T o 2 / 2 M 2 + m 2 + 27tfra cos x 
X = —7= (jSto) exp - /3mc - /3rau —- 

V* L U« + m) 2 

. _ , m + M cos x .1 r / 20Mm , . ,\ 

+ 2fimcu r-^— cos $ \-h ( ■=7-^ cu sin x sin ^ I 

ilf + m J \ilf + ra / 

Passing now to the right hand side of (Ilia) we may replace in the 
first place dC by dy, because of (111b). c and C are then replaced by 
the expressions 

M M 

c = c - tf- i — r + 



M + m * M + m l 

M _ M , 

C M + m r M + m y ' 

With the angles defined in Fig. 15, we thus get for Y 
F = I ^— J (?— 1 exp [— fiMc - j3m(c — u ) ] 

• if [J y 2 d? sin tf0 d8 d<f> 

exp — /3.I/7 2 + 0Mcy (cos ^ cos 6 -+- sin ^ sin cos 6) 

— 2/3 -r-p— 1/7 (cos — cos cos x — sin sin x cos <f>) 

M + m J 
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Integrations over 8 and </> again go with (112). Beofore writing down the 
result we shall pass over to a cylindrical coordinate system defined by 

7ll = 7 cos 6 y± = 7 sin 6 

7 2 dy sin 9 dd = y± dy± dy± 

The result of the first two integrations then reads 

Y = - /3 3 (Mm) 3 ' 2 exp [- fiMt - (Jw(c - u') 2 ] 

IT 

/+°° / m yf \~| 

d-yii exp - fiMy\\ + 2j8M7 M ( c cos ip - M , m ( l ~ cos x) 1 

f 00 / M.7YI i \ 

7 j. (?7x exp [/3M7i]-7 (2/3Mc7x sin ^)-7 I 2/3 ^-p — m 7 x sin x 1 

The first of these two integrals is elementary, the other is Weber's second 
exponential integral 26 which equals 

[ exp {- V Y)h{at)Ubt)t dt = ± exp (£jr) Io (^) ( ll4 > 

This yields for Y exactly the expression (113). The identity (111) is 
thus proved, and with it the convolution theorem. 

The theorem just proved reduces the velocity distribution for arbitrary 
field and temperature to two components, one containing the field, but 
not the temperature, the other the temperature but not the field. In each 
of these components, in turn, the variable parameter scales out; thus the 
general distribution reduces to two basic ones one of which is the 
Maxwellian one: the other is worked out partially in the calculations of 
the Sections IIC and IID. The special case of heavy ion mass has been 
published independently by Kihara 27 without any apparent knowledge 
of this theorem which was available in the literature without complete 
proof. 1 Kihara's form of the theorem is that heavy ions in a light gas 
have an off-set Maxwellian distribution, with the gas temperature as 
parameter if the mean free time condition is obeyed for their collisions. 
Such a function is indeed the convolution of a Maxwellian distribution 
and the 5-function discussed in the Sections IIB and IIC. 

The general distribution function resulting from (107) cannot be 
written down explicitly because this goal was never achieved for h(c). 
However we do find a result which is almost a full substitute for this, 

20 Watson, G. N., A Treatise on the Theory of Bessel Functions. Cambridge 
University Press, Section 13.31, 1922. 

« Kihara, Taro, Rev. Mod. Phys., 24, p. 45, 1952. 
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namely that all averages of products of integer powers of the Cartesian 
velocity components, which were shown to be computable in the high 
field case, can be computed for the intermediate and low field range as 
well. The calculation proceeds as follows. Suppose we wish to compute 
the velocity average 

(c'Xc?) = f c/Vc//(c) dc (115) 

for m, n, p integer or zero. We apply the convolution theorem (107) to 
/(c), decompose the three factors into 

C," = [« a + (ex- u x )} m 

Cv n = K+ (c„ - W„)} n 

Cm' = {lU+ (C.- U,)}' 

and expand each of them by the binomial theorem. We find 



«*w-5s£(:) (;)(;) 

J h(u) u/uyW du f m(vK m ~V~'^ P ~ I dv 



(116) 



The second integral is a thermal average, the first a high field average 
computable by the method of Section IIB. Thus the average (116) is a 
finite sum of products of computable averages and is itself computable. 
When formula (116) is applied to the averages (52), (53), (54), (57) 
and (59) very simple results are found because of the symmetry of the 
function m(v). For the drift velocity (c,) we get from (52) 

(r \ _ M + m // l - cosx\ Ml7 v 

™ ~ —W~l \ ar I (U7) 

This is the same formula as (52) which is thus proved to hold inde- 
pendently of the gas temperature. In the energy formulas we find simple 
addition of the thermal and high field values because the middle term in 
(116) drops out by symmetry. Inserting (53), (54), (57) and (59) we find 

A™,A qz.t j_ CM + mf //l - cos x\ 2 /. 1c n 

{mc ) = 3LT+ ___/ v __) dig) 

(M + mf ( — Si " 2 X + 4m(1 - C0S x) , 

(mcl) = kT + r - i — ™ '— (H9) 

, , 2 /3M sin x + 4m(l — cos x)\/l — cos xN 



\ & A" 



ar 
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(mcl) — m(c z ) 2 = 

,.. . N2 /M sin 2 x + 2m(l - cos x)\ 

(M + m) {— ar i (120) 

,, IZH sin 2 x + 4m(l - cos x)\/l _ cos A* 
3/ \ « A ar / 

> 3 /sin 2 x\ 

, * (tf + m) \-iTZ ( 121 > 

<mc> w -r /3 ^ sin 2 + 4m(1 _ CQS y)yi _ C0B ^\« 



3/ 



\ ar A at / 



The interpretation of these formulas is implicit in the discussion of the 
high field formulas given earlier. In particular the combination of the 
equations (55) and (116) can be given the elegant form 

m(c 2 ) = M{(?) + m{c? + M(c z f (122) 

It states that the energy of an ion is obtained by adding the energy of a 
gas molecule, the energy visible in the drift motion and a storage term 
which is M/m times the energy in the drift motion; this term becomes 
important for electrons in a gas. A low field approximation to this for- 
mula (in which the second term on the right may be neglected) has been 
published in the article of Kihara. 

IIIB. RESULTS FOR THE POLARIZATION FORCE AND THE ISOTROPIC 

"maxwellian" MODEL 

The polarization force between ions and molecules which predominates 
over other forces at sufficiently low temperature satisfies the mean free 
time requirement of the preceding section. It follows that the complete 
theory given for those conditions applies to this force. The magnitude of 
force was given in (4). Its potential equals 

V = - — (123) 

2 p< 

Classical theory is usually applicable to the scattering by the potential 
(123) because angular momentum quantum numbers run as high as 30 
or 50 in normal situations. 29 This classical type theory, first developed by 
Langevin, 30 follows standard elementary methods for computing the 

28 Reference 27, formula 5.12. 

29 Holstein, Theodore, private communication, see also Reference 11. 
"Langevin, Ann. de Chim. et de Phys., 5, p. 245, 1905. 
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angle of deflection x due to a potential of the type (123). The result is 



-*- 2 Cn 



Ul du 

.j , e 2 P(M + m) 4 \ I/2 (124) 



Here b is the "impact parameter", and Wi is the lower of the two positive 
roots of the polynomial in the denominator; if the polynomial has no 
real root, the integration goes from to » . The question whether the 
denominator has a real root or not is tied up with the nature of the orbit. 
If b is sufficiently large a root exists and the orbit looks like a hyperbola, 
Fig. 2(a) ; for small b no root exists and the two particles are "sucked" 
towards each other in a spiralling orbit as shown in Fig. 2(b). The two 
regimes are separated by a limiting orbit in which the particles spiral 
asymptotically into a circular orbit. This limiting orbit is found by 
setting the discriminant of the square root in (124) equal to 0. We find 

6,im = Mm* (125) 

From this value of 6 )im a cross section and a mean free time t, 
for spiralling collisions can be derived. We find 

1 / Mm V" ri9fi s 

T ' = 2^N {{M + m)p) (126) 

This is indeed a constant mean free time as stated, the speed of encounter 
7 having dropped out. 

1/t is the dimensional quantity entering into the averages (y(x)A) 
which occur in the Sections IIB and IIIA. In working them out in detail 
as was done by Hass6 9 one has to take into account hyperbolic collisions 
also; for them a t cannot be defined or comes out to be zero in the mean. 
This is due to small angle deflections which are infinitely probable. 
However, any quantity <p(x)A to be averaged in a physical problem con- 
tains a y?(x) which vanishes for such impacts. Hence finite averages 
result which do not give overdue weight to these types of collisions. 
Following Hass6 9 , we do this in the following way for the present case. 
We write (124) in the form 

r Vl dv 

1 " v + w 
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Here v equals bu and the parameter /3 equals b/b lim . It is the parameter 
/3 introduced by Hasse\ Now by the definition of t we have 

M^\ = N J y<p(x) d*(y,x) 

= N [ y<p(x)b db f de 
Jo J o 

m rNyb\ la [ tp(x) d((3 2 ) 

Jo 

From (125) and (126) the factor in front of the integral just equals 1/t, ; 
the integral on the other hand is a computable pure number independent 
of 7 which is obtained by inserting into it the relationship (127) between 
X and &. Hence we may write 

Mx)\ = i_ r M dW (128 ) 

\ T I T s Jo 

The three equations (126), (127) and (128) completely define the nature 
of the averages appearing in previous sections. The integral (128) has 
to be computed by numerical methods. It is seen in the course of the 
evaluations that it naturally decomposes into two parts. The part for 
which /3 varies from to 1 deals with spiralling collisions and exists for 
any p(x). For between 1 and » we get the contribution of the hyper- 
bolic collisions to the average. This part is only finite if ^(x) vanishes for 
small angle deflections. 

The averages (52), (53), (54), (57) and (59), as well as (117) to (121) 
contain numerous averages of the form (128) all of which satisfy the 
predicted condition <p(0) = 0. They are obtained by linear combination 
of two basic types: ((1 - cos x)A> and (sin 2 x/t). The first average is 
given in Hasse\ 9 Separating the parts due to spiralling and hyperbolic 
collisions we find 

f (1 - cos x) d(p 2 ) = 0.8979 
Jo 

f° (1 - cos x) d{^) = 0.2073 

This combines to give 

( l ~ C0S X ) = 1-1.1052 (129) 
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The analogous result for sin 2 x was obtained by the computing group 
of Bell Telephone Laboratories 

I sin 2 x d(tf) = 0.511 

/CO 
sin 2 x d{f) = 0.261 

(™L*) = 1.0.772 (130) 



which gives 



We may now rewrite the major results of Section IIIA for the polariza- 
tion force. From equation (117) we get 

, v 0.9048 . /I " 1 E 

<Cz> = -%r V m + « *v? (131) 

This formula may be found in the literature. 31 What is new about (131) 
is the realization that it is exact at high as well as low electric field. 

The formula for the total energy needs no discussion for a special 
model; it does not involve the angular distribution law when written in 
the form (122). Thus we would obtain, for instance, for an ion travelling 
in the parent gas that its total energy is obtained by doubling its ap- 
parent energy observable in the drift and adding to this the thermal 
energy %kT. 

For the partition of the high field component of the energy in the 
three coordinate directions we have two formulas, formula (58) parti- 
tions the entire field contribution of the kinetic energy, formula (60) only 
its random component. The first formula gives 

e x :e v :e t = M:M:(M -f 6.73m) (132) 

Formula (60) gives 

e x :e v :e* = (M + m) : (M + m) : (M + 3.72m) (133) 

It is convenient to apply the general formulas also to the case of con- 
stant mean free time, coupled with the assumption of isotropic scattering. 
This combination of assumptions represents, strictly speaking, an im- 

81 The formula ia equation (3), p. 39 of Reference 2, in the limit X = 0; or also the 
last unnumbered equation on p. 919 of Reference 6. 
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possibility; for we know of no mechanical force which realizes this ar- 
rangement. This model was already taken as the basis of the Monte Carlo 
calculation in Section IID. It will be seen now that it has a wider sig- 
nificance than one might anticipate. The necessary angular averages are 

(1 - cos x > = 1 (134) 

(sin 2 X > = % (135) 

This yields for (117) 

(O - M or d 36 ) 

As usual, the formula for the energy does not involve the law of scattering 
if written in the form (122). If we choose the form (118) instead we get 
in agreement with (79) 

2 ._ (ilf+m) 3 2 2 fio 7 \ 

(mc ) = SkT H irF z — a r Wit) 

The partition formula (58) becomes 

6,16,16, = M:M:(M + 6m) (138) 

the partition formula (60) which counts random energy only becomes 

e x :ey:e* = (M + m):(M + m):(M + 4m) (139) 

Comparison of these expressions with the ones for the polarization 
force shows that the difference between it and the isotropic model is 
remarkably small from a kinetic standpoint. We may see this by com- 
paring (132) and (138) or (133) and (139). For the other formulas, we 
may compare more specifically the polarization results with an isotropic 
case having its mean free time t given by 

r = 0.9048 r. (140) 

Equation (136) becomes then identical with (131) and because of (122) 
the same identity persists for the energy formula (137). In the light of 
this we may say that it is very nearly correct to state that scattering is 
isotropic for the polarization force. This qualitatively correct fact was 
repeatedly made use of in the preceding sections of the paper. The 
reason for it is chiefly the predominant effect of spiralling collisions. 
Indeed, equation (140) shows that a modification of t, by only 10 per 
cent takes into account the main influence of hyperbolic collisions. 
From the discussion in Section IA it may be seen that the results 
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obtained for the polarization force have a potentially wide field of ap- 
plication when measurements of ion drift are extended to low tempera- 
ture. In the meantime, the results apply occasionally at room 
temperature, whenever we deal with a small ion and are not bothered 
by special scattering mechanisms having large cross section. An example 
of this are the molecular noble gas ions in the parent gas whose drift 
velocities were measured by Hornbeck 16, ' and Varney. Table IV shows 
the measured mobility at standard gas density measured for these ions, 
in comparison with a value obtained from equation (131). The field 
range from which the observed mobility was obtained is intermediate. 
There is not only good numerical agreement, but the experiments follow 
the theory also in that there is little variation of the observed value 



Table IV 

Mobilities at standard density of the noble gas molecular ions. Comparison of 
the experiment with a formula based on the polarization force only. 



Gas 


cm 1 


cm 1 


Mob '- Volt sec 


M€a,B Volt sec 


He 

Ne 

A 

Kr 

Xe 


18 
6.5 
1.9 
1.2 
0.7 


18.2 
6.21 
2.09 
1.18 
0.74 



with the field. The discrepancy between the two columns can be used to 
determine a hard collision cross section which is to be superimposed on 
the polarization force, as is suggested in the so-called Langevin model. 

IIIC. VELOCITY DISTRIBUTION FUNCTION FOR ELECTRONS 

We have almost exhausted the results achieved for intermediate 
field conditions. For the sake of completeness I shall mention shortly 
the intermediate field distribution function for electrons whose derivation 
we owe to the ingenuity of Davydov. 

The derivation does not differ in principle from the one presented in 
Section IIC for the electrons in the high field case. The distribution 
function is first expanded in spherical harmonics. For group theoretical 
reasons the scattering term in the Boltzmann equation is diagonal in 

32 Davydov, B., Phys. Zeits. Sowjetunion., 8, p. 59, 1935. See also Reference 4, 
pp. 349-350. 
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such a decomposition even in the presence of molecular agitation. Thence 
a generalized form of (47) may be derived containing essentially the 
.same terms. Finally, all but the first two spherical harmonics are dropped 
and two equations analogous to (63) and (64) are obtained. In fact, it 
is found that equation (63) is maintained entirely. An extremely compli- 
cated reasoning is required, on the other hand, to find the generalization 
of (64). The result is 

*»-{¥S+»3«»)fcS^ (141) 

Combining (63) and (141) we find 

1 3kT\ df Q . n m 



/I - cos x\ + W ) fc + 3 M Cfo ~ ° 
\ ar(c) I J 



and hence 



/o(c) = exp 



— m I 
Jo 



c dc 



k" 



/ l - cos x \ 2 
\ ar{c) I 



4- kT 



(142) 



This is the so-called Davydov distribution which is a generalization 
containing within itself the Maxwellian distribution as well as the high 
field distribution (65). 

The mean energy and the drift velocity of electrons may be calculated 
from (63) and (142). They are obtainable from the literature and will 
not be discussed here any further. Equipartition of the energy exists at 
all field conditions. 

Part IV — Diffusive Motion of Ions 

IVA. DIFFUSION FOR MEAN FREE TIME MODELS 

It was proved in Section IC that if there are spatial inequalities in the 
distribution of the charge carriers then a smoothing out process sets in 
which can be described as diffusion. This derivation of principle can be 
supplemented for "Maxwellian" molecules by an explicit computation 
of the two components of the tensor (24), that is an evaluation of the 
integral (23). We shall do this by following the method of Maxwell 19 
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rather than by generalizing the formal procedure of the Sections IIA, 
IIB and IIIA. Such a generalization would no doubt be possible, but 
would increase unduly the bulk of this paper. We shall operate therefore 
directly on equation (20). To get out the integral (23) we multiply the 
equation vectorially with c and integrate over dc. This operation makes 
the first term vanish completely. This is obvious from symmetry for the 
components c x and c„ of the multiplier c. For c, we have 



•/*• 



dc 



An integration by parts brings this in the form (18) and thus makes 
it equal to zero. 

Temporarily, we may break the integral term of (20) into two parts, 
using some artificial procedure to eliminate small angle collisions. The 
first half of the integral term reads then simply 



- J g(c)c dc 



This is already the desired average (23). On the second half we use the 
identity (108} to give it the form 

-i // W)0(c')cll(x) dO, dC' dc 

We now use (7) to replace c by the expression 

m / . M _, . M 

c = c ^ C 4 y 

M + m M + m T M + m T 

Only Y i s affected by the integration over dfi 7 which we take up first. 
Using y as the axis of a polar coordinate system we may write 

Y = Yll + Y-<- 
For every value of x, Yi I has the fixed value y cos x- On the other hand 
the average of y-j- vanishes through integration over all azimuths. Thence 
we may write 

CII(X) dQy = -4— C ' + j^r— C + tJL- ( C ' - C') (COS X> 

= m + M (cos x) / . Af(l — cos x) «/ 
M + m C M + m 

We now multiply with M(C')g(c') and integrate over dc' dC'. The inte- 
gration of the term containing C obviously vanishes for two independent 
reasons. The integration of the term in c', finally, yields again the average 



4tt J 
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(23) . Combining the two pieces, we find 

M ( l -*<**) f g( C ) C dC 

M + m\ r / J " w 
In this expression, the artificial exclusion of small angle scattering is no 
longer necessary and can be dropped. Completing the integrating of 
equation (20) we see that the right hand side gives averages over the 
unperturbed velocity distribution /(c). Combining pieces, using (23) 
and indices 1, 2, 3 for the x, y and z components we get 



a 

tit, t) = -n(r, t) £ k. 



M + m 



M 



/l - cos x\ 



Uac) - (ci)(c r )} 



(143) 



\ — ; — / 

According to (16) and (24), the square bracket in (143) is the diffusion 
tensor. It has two distinct components which equal respectively 

M + m (c?) - (c z ) 2 



D„ = 



Dx - 



M 



/l — cos x\ 

\ — ; / 



M + m (ex 2 ) 



M /l — cos x\ 

V 



(144) 



(145) 



T / 

The velocity averages entering are (120), and (121), that is the directional 
components of the random part of the energy. Substituting we get finally 

(M + m)kT 



Dm = 



Mm 



+ a 2 



/l - cos x\ 

\ r / 

... x ,3 /M sin 2 x + 2m( l - cos x) 2 \ 
(M + m)\ T l 

/ZM sin 2 x + 4m(l - cos x)\A - co s x\ 3 



(146) 



M 2 m 



V 



A" 



7 



Dx = 



(M + m)kT 
, /l — cos x\ 



+ a 2 



- N 4 /sin 2 x\ 

(M + m) \— ^— / 



(147) 



, ro /3Af sin 2 x + 4m(l - cos x )\/ l - cos x \ 3 
M*m \ ; /\ ~ / 



The diffusion coefficients have the simple property that they are ob- 
tained by adding the low field and the high field limiting expressions. 
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This is a consequence of the limited form of the convolution theorem 
proved in Section IIIA; it probably implies also that the theorem can 
be extended in some form to include the case of diffusion. 

It has been mentioned in the Section ID that the Nernst-Townsend 
relation (30) applies only to ions moving in a low field. We are now in a 
position to examine possible extensions of it to general fields. Equations 
(144), (145) and (117) suggest the form 

Dn _ 2 X mean random energy along n . . 

mobility e 

where n stands for one of the principal directions of the diffusion tensor. 
This formula contains equation (30) as a specialization to the low 
field case. 

Formula (148) is one of the formulas obtained in this study of ion 
motion in which model parameters do not appear. It is valid (a) for all 
interactions at low field and (b) for the mean free time case at all fields. 
It also holds dimensionally at high field for models obeying (25); this 
may be seen from (26a) and (28a). It appears a reasonable conjecture 
that (148) is approximately true for any law of interaction; the question 
will be taken up again in the next section. 

Let us, in conclusion, write down the formulas resulting from (146) 
and (147) for the two special mean free time models studied in detail 
in Section IIIB: the polarization force and the isotropic model. The 
necessary averages are (129), (130), (134) and (135). They yield for the 
polarization force 



(149) 



Z>„ = M J" m 0.905iv feT 

Mm 

1 (M + m)\M + 3.72m) t( , 

+ 3 MHn{M + 1.908m) fl ( °- 905t ' ) 

and for the case of isotropic scattering 

M + m , T 1 (M + m)\M + 4m) 23 , 1K1 v 

D " = -Mm~ TkT + 3 Mhn{M + 2m) " T (151) 

M + m _ hrp , 1 (M + m) 4 _, , 
D± ~ ~Mm- TkT + Z M>m(M + 2m) aT (152) 

Just as in the earlier study the results for the two models do not differ 
appreciably. 
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IVB. LONGITUDINAL DIFFUSION FOR THE HARD SPHERE MODEL 

Whenever the mean free time condition for collisions is not fulfilled, 
then the computation of diffusion coefficients requires a procedure analo- 
gous to that of Section HE. Since this entails some numerical work the 
calculation was only carried out for a case which was thought to be of 
experimental interest, namely for longitudinal diffusion of ions in the 
parent gas. In other words, we are extending the numerical computation 
at the end of Section HE to include longitudinal diffusion. The computa- 
tion to provide us with the undetermined constant of equation (28b) 
for the special case when m and M are equal; it also offers, incidentally, 
a good test case for applying the method of Section HE outside the 
area for which it was designed originally. 

Since the equation is only to be solved in the high field case we may 
apply to (20) the reduction method of Section HA. If we introduce also 
the specialization warranted by the hard sphere model and unit mass 
ratio then, in analogy to equation (40), we get the following starting 
equation 

dtf(w) . ( -v 1 C w'* sin x dx f * „/ _/\j„ 

1ST + wg(w) _ 5 J. — ^— I. ° {w ) * (153) 

= — \k{w z — (w z )}h(w) 

Here the dimensionless variable w defined by (85) has been employed 
instead of c. 

Equation (153) is the fundamental equation of our problem; it is an 
inhomogeneous version of equation (40). We solve the equation in the 
same way as we did previously, namely by decomposing g(vt) into spheri- 
cal harmonics and forming moments. In other words we follow step by 
step the procedure of Section HA, the only difference being the presence 
of an inhomogeneous term. We shall not enumerate all these steps again. 
We shall only note in passing the inhomogeneous form of (47) which is 

- / ^r g>(w')P> (cos k) sin x dx - wgXw) 
2 Jo w 3 

dg^ v-l„ _ i(w) 



2v — 1 1 dw w 



2v + 3 ( dw w j 

= M L - , wh^i(w) + "* \ wh, + i(w) - (w,)h(w) 
\_2v — 1 2v + 3 J 
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Having introduced moments in the manner described earlier we arrive 
at the inhomogeneous version of (49) or (86) 

v(8 + v + l){s - 1, v - 1} + („ + l)(s - v){s - 1, v .+ 1) 

-(2, + 1)(1 - (/..,»{s + 1, p) = -(2, + 1)<1, 1)< 8 , „) (154) 

+ !•<• + 1, F - 1> + (F + 1)<S + 1, » + 1) 

Here the curly brackets [s, v) are normalized moments over g(w) defined 
as follows 

{«, v) - £- J giyr)w'P r (cos tf) dw (155) 

The equations (154) show that the quantities {s, ?} are numbers, the 
variable density gradient nk having been eliminated by the definition 
(155). The system does permit that arbitrary amounts of the pointed 
averages be added to the curly ones. This indeterminacy is removed by 
the supplementary condition (18) which, in the present notation reads 

[0, 0} =0 (156a) 

The connectivity of the equation system (154) is the same as that of 
(86). Hence it will have the same properties as that earlier system. We 
may, therefore, reduce it in the manner followed previously and get 
inhomogeneous versions of the equations (87), (88) and (89). They read 

{4,0} = -jj (4, 1) + | <1, 1)(3, 0). (157a) 

-y<2,0)-^(2,0>+10<l,l> 2 
112(3, 0} - 3{7, 0} = 4(7, 1) - 4(1, 1)(6, 0> 
+ f (5,0)4-1^(5,2) 

- l f (1, 1X4, 1) 



1344 /q i\_i_ 1344 /q o\ (158a) 



25 * ' ' ' 25 

+ ^ (1, 1>(2, 0> - ^ (1, 1>(2, 2) 



250 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1953 

54{2 , 0} - f {o,o} + £ {io,o } - -^ <io, i) + "a 1><9,0) 

-g<8,0>.g<8,2) + H (1 ,l )(7 ,l) 

+ 6(6) 1) - 2 1 (6, 3) - ^ (1, 1X5, 0) + 7(1, 1X5, 2> (159a) 
+ ^ (4)0) + 15_ 8 ( 4 ,2>-f(4,4> 

-1|? (1,1X3,1) + Z?<1,1X3,3> 

The pointed averages over the distribution h(w) may be found in Table 
II. Substituting them we get 

{0, 0} = (156b) 

{4,0} = -10.494 (157b) 

112{3, 0} - 3{7, 0} - 647.8 (158b) 

54{2, 0} - ||5 {6, 0} + Jl {10, 0} = -566.4 (159b) 

The form (90) that was assumed for fo(w) will again be taken for g (w) 
with new undetermined coefficients p, q, r, s and a factor k\ evident from 
(153) or (155): 

g (w) = k\[ V Ei{\v?) + qlU{W) + *** + *U?Kx{W)\ (160) 

This is a rather poor assumption because the form (90) was adopted for 
ho(w) after an extensive study of the properties of the distribution 
function /i(w). For ^r(w) we know little beyond the fact that it is some 
kind of distorted p-type function. The g (w) derived from this is not 
likely to resemble h (w) very closely. Thus the choice (160) is mainly 
based on ignorance and convenience; this explains the slower convergence 
observed here than in (91), (93) and (95). To start with, the zero order 
is completely lost because (156) yields a zero coefficient. We find in first 
order, using (156) and (157) 

p (1) = 4.8842 a) _ (i) _ n 

g (1) = -4.2689 r " S 
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in second order, using (156), (157) and (158) 

p (2) = -10.542 

q m = +14.993 s (2) = 

r (2) = -4.408 
in third order, using (156), (157), (158) and (159) 

p (3) = -0.8710 

q (3) = 4-1.1754 

r m = +1.0140 

s (3) = -0.5809 

The longitudinal diffusion coefficient results from these numbers by the 
use of (23), (24) and (160). With the notation (155) the formula becomes 

D u = -«W»|l f l) (161) 

The formula (154) yielding {1, 1} from g (w) is s = 2, v = 

[1, 1} = |{3, 0} + i(3, 1) - *<1, 1><2, 0) (162a) 

or numerically from the Table II 

{1,1} = i{3, 0} +0.6433 (162b) 

The result is 

{1, 1} (1) = -0.3695 (163a) 

{1, 1} (2) = -0.2075 (163b) 

|1,1} (3) - -0.2198 (163c) 

The numbers do not extrapolate too reliably but one would guess that 

{1, 1} = -0.22 

is essentially correct. Hence we have 

2>l, = 0.22a 1/2 X 3/2 (164) 

In order to gain an appreciation of the value obtained it is worthwhile 
to compare it with the value that would have been predicted from the 
generalized Nernst-Townsend relation (148). The mobility concept is 
ambiguous for all but the cases discussed then. It would seem that the 
appropriate concept here is the differential mobility because comparison 
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is made between a small density gradient and a small change in the 
applied field. Thus we would interpret (148) to mean 

da 
which, with (85), (92) and (96), becomes 

2>„ « 0.26a 1/a X 3/2 (165b) 

The error of formula (165) is thus 18 per cent, when compared to (164). 

Part V — Concluding Observations 

The present article is supposed to contain the essentials of a kinetic 
theory of charged particles moving through a gas in the presence of an 
intermediate or high electric field. An effort was made to make the theory 
general, yet many irksome restrictions will become apparent to those 
who will try to apply it to their particular problem. Especially those who 
have in mind application to electrons will find the article unsatisfactory. 
It is true that many sections leave the masses variable; however, the 
assumption of elastic collisions, which is made throughout, is almost 
fatal to all but the most elementary applications. Thus most of the 
material is slanted for ions. Within this domain, numerous awkward 
restrictions are still found here. The most important ones are pre- 
sumably the restriction to D.C. conditions, the assumption of "low" 
ion density, and the omission of all magnetic effects. It is my general 
impression, which I gained from the convolution theorem Section IIIA 
and which is confirmed by a recent publication 27 that much can be done 
to remove these three restrictions provided the mean free time assump- 
tion is made for collisions. To many the adoption of the mean free time 
condition will in itself appear an awkward restriction. In a rigorous sense 
this is true, and calculations are made in this article for the more ap- 
propriate hard sphere model when quantitative comparison with experi- 
ment is contemplated (equations (100) and (164)). Indications are even 
given for a treatment which dispenses altogether with the use of models 
(equation (103)). However, for rapid advance and easy handling, the 
mean free time assumption does appear essential. It is therefore im- 
portant to point out that in a wider semiquantitative sense, the use of 
this model is no barrier to application. In other words, there is in the 
mean free time formulas information which suggests a wider validity. 
This is particularly true for equations which do not contain model 
parameters, such as (55), (56), (122) and (148). Even formulas which 
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do contain the mean free time yield to judicious treatment. For example, 
we have the hard sphere formula (100) for the drift velocity of an ion. 
This formula happens to be limited to the high field case and mass ratio 
unity. On the other hand we have formula (136) which holds for all fields 
and all mass ratios, but assumes constant mean free time. We now adopt 
this formula as a general guess for the hard sphere model, interpreting 
r as previously as the mean free time between collisions; this quantity 
is now no longer a constant, but should be taken as 

T - VW+ (C>) (166) 

The denominator is the root mean square relative velocity which is 
familiar from other applications. The interpretation (166) yields a tenta- 
tive formula for the drift for all mass ratios and for all fields. Specializing 
to the high field case, we may neglect (C 2 ) in (166) and then substitute 
for (c 2 ) from (55). This yields the high field formula 

(^ ( ^;»" v " tf d67) 

This is indeed a very successful formula. For ions in the parent gas it 
differs from (100) by only 4 per cent. For electrons it checks the result 
of Druyvesteyn to within 12 per cent. Finally, for heavy ions in a light 
gas, we find exact agreement with equation (71). As a second specializa- 
tion we may apply (166) to the low field case. We must then set 

« + <*> -**(!+£ 

and get from (136) and (166) 

1/1 1V /2 eEX , x 

All dimensional factors in this formula are correct. Numerically (168) 
is somewhat inferior to (167); for the factor differs from the correct one 33 
by 20 per cent. Nevertheless, the combination of (136) and (166) gives 
results which are semiquantitatively correct in all relevant limiting cases. 
This makes it a reliable interpolation formula for intermediate field con- 
ditions; for this case (c) would have to be substituted from (122) and 
the resultant quadratic equation solved for (c 2 ). 

From the examples given we may conclude that the mean free time 
formulas contain in essence information applicable to other types of 
elastic scattering. 

35 See Reference 2, page 40, second equation. 
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